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CHAPTER  1:  THE  FHENCMENOIDGY  AND  TERMDTOIDGY  OF  OQSMTC  BAY  ACCESS 

1.1  PREFACE 

This  chapter  presents  a  summary  of  the  current  understanding  of  the 
phenomenology  of  cosmic  ray  cutoffs.  A  list  of  tens  is  presented  fay  which 
cosmic  ray  cutoffs,  and  access  generally,  may  be  discussed.  The  text, 
which  originally  constituted  an  as  yet  unpublished  paper  "Cosmic  Bay 
Cutoffs  Within  Planetary  Magnetic  Fields  *  An  Examination  of  the 
Phenomenology  and  Terminology*'  prepared  by  the  author  of  this  report  and 
eight  other  researchers  (listed  below) ,  grew  out  of  the  papers  "On  the 
Terminology  of  the  Penumbra"  by  Oooke,  Shea  and  Smart,  1981;  and  "Be- 
evaluation  of  Cosmic  Bay  Cutoff  Terminology"  by  Cooke  et  al.  (the  nine 
researchers) ,  1985. 
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Nucleaires  de  Saclay,  Sad  ay,  France 


l 


1.2  INTRODUCTION 


The  question  of  charged  particle  access  to  points  within  the 
geomagnetic  field  has  been  a  topical  one  throughout  the  entire  history  of 
cosmic  zay  physics.  Results  obtained  from  theoretical  investigations  into 
the  phenomenology  of  charged  particle  access  have  been  instrumental  in 
aiding  the  interpretation  of  a  vide  range  of  experimental  investigations. 
These  studies  have  ranged  from  examination  of  the  aurora,  through  attempts 
to  account  for  the  observed  directional  asymmetries  detected  in  the  primary 
and  secondary  oosmic  ray  fluxes,  cosmic  ray  astronomy,  the  determination  of 
primary-secondary  oosmic  ray  relationships,  to  the  utilization  of  specific 
access  properties  of  the  geomagnetic  field  for  the  determination  of  oosmic 
ray  isotopic  abundances. 

There  is  a  growing  interest  in  the  question  of  access  of  oosmic  rays 
into  the  magnetospheres  of  other  planets  (Jupiter  is  of  particular  interest 
at  this  time) ,  specifically  in  regard  to  the  charged  particle  fluxes  likely 
to  be  encountered  by  spacecraft  in  the  vicinity  of  the  planets,  and  in  the 
interpretation  of  the  detail  of  the  flux  distributions  actually  measured, 
and  the  information  these  yield  on  the  characteristics  of  the 
magnetospheres. 

The  early  workers  who  researched  the  question  of  charged  primary 
oosmic  ray  particle  access  to  positions  within  the  geomagnetic  field 
recognised  that  the  general  equation  of  charged  particle  motion  in  a 
magnetic  field  does  not  have  a  solution  in  closed  farm,  and  developed 
alternative  techniques  for  investigating  the  aooess  problem.  The  initial 
papers  in  this  area  were  published  by  Stormer  (1930) ,  Bouckaert  (1934) , 
Lanaitre  et  al  (1935)*  and  Lamaitre  and  Vallarta  (1936a, b)*.  The 
attentions  of  those  authors  were  directed  towards  the  need  to 
systematically  distinguish  between  directions  of  arrival  along  which  a 
particle  from  interplanetary  space  could  reach  a  given  location  (allowsd 
directions)  and  those  directions  not  accessible  to  such  partidss 
(forbidden  directions) . 

The  approach  of  these  pioneering  investigators  took  advantage  of  ths 
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fact  that  the  magnetic  field  of  the  Earth  could  be  approximated  by  that  of 
a  sisple  dipole.  In  such  an  axially  symmetric  field,  bound  periodic  orbits 
are  especially  significant  in  delimiting  and  characterizing  the  different 
access  regions.  (The  bound  periodic  orbits,  which  are  those  totally 
constrained  to  remain  within  the  magnetic  field,  possess  a  regular 
periodicity.  This  periodicity  is  manifested  as  cyclic  variations  in 
altitude,  latitude,  and  longitude  along  the  trajectory.) 

Experimental  investigations  of  the  access  problem  were  carried  out  by, 
amongst  others,  Malmfors  (1946) ,  Brunberg  (1953) ,  and  Bruriberg  and  Dattner 
(1953) .  These  authors  used  experiments  in  which  the  motion  of  electrons  in 
the  proximity  of  a  charged  terella  was  used  to  model  cosmic  ray  motion  in 
the  geomagnetic  field. 

With  the  advent  of  fast  digital  computers  it  became  practicable  to 
establish  the  extent  of  significant  allowed  and  forbidden  regions  by 
systematic  direct  checking  of  access  of  individual  trajectories  calculated 
using  high  order  simulations  of  the  geomagnetic  field.  The  technique  was 
immediately  applied  to  confutations  which  were  beyond  the  practical 
capability  of  the  early  investigators,  such  as  determination  of  asynptotic 
cones  of  acceptance  for  the  world  wide  network  of  cosmic  ray  stations 
(McCracken  et  al. ,  1962,  1965,  1968)  and  calculation  of  a  world  wide  grid 
of  vertical  cosmic  ray  cutoff  rigidities  (Shea  et  al.,  1968). 

In  these  initial  digital  computer  studies  it  was  unnecessary  to 
correlate  the  position  or  extent  of  the  allowed  and  forbidden  access 
regions  with  the  associated  trajectory  characteristics.  As  a  result,  far 
many  years  there  was  little  effort  to  relate  the  real  field  access  patterns 
to  those  predicted  for  the  sinple  axially  symmetric  field,  or  to  look  into 
the  role,  if  any,  that  bound  periodic  orbits  play  in  defining  the  allowed 
and  forbidden  regions  in  the  reed,  field.  Nevertheless,  the  terms  used  in 
the  dipole  case  to  refer  to  the  access  regions  were  carried  ever  into  the 
numerical  calculations,  although,  as  noted  by  Cooke  (1981)  and  Gall  (1981), 
often  in  an  imprecise  way.  Whereas  the  early  theoretical  workers  viewed 
access  conditions  in  what  may  be  called  a  "direction  picture",  describing 
the  directions  from  which  particles  of  a  specific  rigidity  oculd  or  could 
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not  arrive,  the  later  cxcputer  calculations  have  usually  used  a  "rigidity 
picture"  in  which  accessibility  is  considered  as  a  function  of  particle 
rigidity  in  a  single  arrival  direction.  The  use  of  the  old  geometric 
terns,  appropriate  to  the  direction  picture,  has  caused  considerable 
confusion  in  the  rigidity  picture  studies. 

It  is  believed  that  a  reexamination  at  this  time  of  the 
characteristics  of  the  access  regions  distinguished  by  Sterner  and  by 
Lemaitre  and  Vallarta,  together  with  a  ccnpariscn  between  the  properties  of 
these  analytically  distinguished  regions  and  those  detected  by  the  digital 
oenputer  method,  will  be  useful  for  further  progress  in  the  understanding 
of  the  cosmic  ray  access  problem. 

Additionally,  the  developing  interest  in  examining  the  characteristics 
of  the  cosmic  ray  fluxes  within  other  planetary  magnetospheres  has  spurred 
interest  in  transposing  the  new  techniques  and  insights  across  from  the 
terrestrial  work. 

This  chapter  reviews  the  past  and  present  work,  and  then  presents  a 
set  of  definitions  which,  it  is  hoped,  will  meet  present  needs  in  relation 
to  the  consideration  of  access  to  locations  within  the  magnetic  field  of 
any  planet. 


1.3;  HISTORICAL  PERSPECTIVE 
1.3.1  Magnetic  Field  Alone 

Stormer  (1930) ,  later  summarized  by  Stormer  (1955) ,  considered  the 
matter  of  aooess  from  the  point  of  view  of  fully  forbidden,  rather  than 
fully  allowed,  directions  of  arrival.  In  addition  to  establishing  soma  of 
the  important  general  properties  of  charged  particle  trajectories  in  the 
geomagnetic  field,  he  showed  analytically  that  at  ary  point  in  a  magnetic 
dipole  field  there  exists  a  right  circular  ocnical  shall  of  directions  of 
arrival,  within  which  all  aooess  is  forbidden  to  particles  of  a  specified 
magnetic  rigidity  approaching  from  infinity. 


The  opening  angle  of  this  acne  is  a  function  of  particle  rigidity,  and 
its  axis  is  perpendicular  to  the  dipole  axis  and  also  to  the  radius  vector 
from  the  dipole  oenter  to  the  point  concerned.  In  terrestrial  terms  the 
axis  of  the  cone  is  horizontal  at  this  point  and  is  directed  east-west. 

The  existence  of  this  ocne,  which  is  often  called  the  Sterner  acne  and 
which  encloses  forbidden  directions  of  arrival,  is  purely  a  property  of  an 
axially  symmetric  field,  and  is  independent  of  the  presence  or  absence  of  a 
solid  body,  such  as  the  Earth,  which  might  envelop  the  source  of  the  field. 

For  positively  charged  cosmic  rays,  all  directions  to  the  east  of  the 
Stormer  ocne  are  populated  only  by  particles  in  trend  periodic  orbits.  As 
particle  rigidity  increases,  the  opening  angle  of  the  cone  decreases.  Par 
any  given  direction,  a  particular  rigidity  value  exists  for  which  the 
Stormer  cone  lies  in  that  direction.  At  this  rigidity  value,  and  at  all 
lower  values,  all  access  from  infinity  in  the  given  direction  is  forbidden. 
In  the  rigidity  picture  this  limiting,  or  "cutoff",  rigidity  value  is 
termed  the  "Stormer  cutoff  rigidity" ,  or  simply  the  "Stormer  cutoff". 

It  is  essential  to  note  that  there  has  been  a  change  in  the  usage  of 
the  term  "Stormer  ocne"  since  it  was  first  introduced.  Early  authors  (e.g. 
Lemaitre  and  Vallarta,  1936b;  Hutner,  1939a)  used  the  term  to  identify  the 
boundary  between  fully  inaccessible  and  partly  accessible  solid  angles  of 
particle  arrival  an  the  direction  sphere.  The  visage  changed  fran  this 
definition  of  a  she!  1-1  Dee  boundary  to  one  incorporating  the  solid  angle 
containing  all  arrival  directions  which  are  not  expressly  forbidden  by 
Stormer  theory  (Alpher,  1950;  Kasper,  1959;  Vallarta,  1961,  1971).  With 
the  advent  of  digital  computer-based  calculations  the  usage  became  reversed 
to  represent  the  solid  angle  which  is  fully  forbidden  (e.g.  Schwartz, 

1959) .  Great  care  must  therefore  be  taken  when  consulting  the  literature 
cn  this  matter.  To  alleviate  further  confusion  the  term  "Stormer  ocne" 
will  be  used  throughout  this  report  in  its  original,  shell-like,  sense,  and 
the  term  "forbidden  ocne"  (section  1.5)  will  be  used  to  describe  the  solid 
angle  within  which  all  access  is  forbidden  to  charged  particles  of  a 
specified  rigidity  and  sign.  In  figure  1.1  the  forbidden  ocne  far 
positively  charged  particles  is  the  solid  angle  region  which  is  epen 
towards  the  east.  Stormer  theory,  pertaining  as  it  does  only  to  the  natter 
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FIGURE  1.1 

Notional  spatial  relations  of  the  allowed  cone,  sain  cone,  penumbra. 
Sterner  acne,  and  forbidden  acne  for  positively  charged  cosmic  ray  nuclei 
with  an  arbitrary  rigidity  value  at  an  arbitrary  location  in  a  magnetic 
dipole  fiaM. 
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of  forbidden  access,  makes  no  predictions  about  the  inportant  question  of 
accessibility  in  directions  outside  the  forbidden  cone.  The  essential 
point  is  that  on  one  side  of  the  Stonier  cone  charged  particle  access  is 
forbidden  for  all  directions,  whereas  on  the  other  side  of  the  cone  charged 
particle  access  may  occur  in  certain  directions. 


Although  Sterner 's  original  work  was  done  for  a  dipole  field,  the 
theory  applies  equally  well  to  any  other  axially  symmetric  field,  and  other 
investigators  have  locked  into  the  characteristics  of  the  forbidden  cane 
for  more  cccplex  representation  of  the  Earth's  field.  In  each  case  an 
additional  ccnplexity  was  introduced  into  the  axially  symmetric  field  in 
order  to  represent  the  field  in  what  was  perceived  to  be  a  more  realistic 
way  (e.g. :  a  dipole  plus  ring  current  -  Trieman,  1953;  dipole  plus 
quadripole  -  Queriby  and  Webber,  1959;  and  dipole  plus  uniform  field  - 
Obayashi  and  Hakura,  1960) .  The  field  of  Jupiter  has  also  proven  to  be 
essentially  dipolar  in  form  (although  the  higher  order  terms  are  found  to 
be  basically  dipolar  in  form  (although  the  higher  order  terms  are  found  to 
be  somewhat  higher  relative  to  those  in  the  geomagnetic  field) ,  and  so  the 
Sterner  theory  can  be  applied  to  the  study  of  the  cutoffs  for  Jupiter  (see 
for  example,  Cooke,  1974) . 


1.3.2  Effect  of  the  Solid  Planet 

Stormer's  work  was  of  profound  importance  for  determining  the 
filtering  effect  of  a  dipole  field  on  the  incoming  cosmic  radiation,  but  it 
did  not  take  into  account  an  effect  which  turned  out  to  be  extremely 
Important  for  cosmic  ray  research  carried  out  close  to  the  surface  of  the 
earth:  the  presence  of  the  solid  earth  itself,  together  with  its  atmosphere 
(or  more  generally,  the  presence  of  the  solid  planet  within  its 
magnetosphere) .  At  rigidities  marginally  above  the  Stormar  cutoff,  the 
trajectories  followed  by  incoming  charged  particles  in  a  dipolar  (or  almost 
dipolar)  field  like  that  of  the  earth  are  often  very  oerplex,  going  throuepi 
a  number  of  altitude  oscillations.  In  the  case  of  the  earth,  if  the 
particle  dives  deep  into  the  atmosphere  during  one  of  these  oscillations, 
it  is  lost  by  interaction  with  atmospheric  nuclei  and  cannot  then  be 
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observed  further  along  its  theoretical  track.  (This  track,  which  we  refer 
to  later  in  this  report  as  a  'Virtual  trajectory",  is  the  track  which  the 
particle  would  have  followed  had  the  planet  not  been  present  and  had  the 
trajectory  been  control  led  solely  by  the  magnetic  field.)  In  practice  the 
presenoe  of  the  earth  excludes  fran  observation  most  of  the  particles  with 
arrival  directions  immediately  to  the  west  of  the  Stormer  cone.  Therefore, 
to  successfully  model  the  cutoff  distribution  pattern  which  actually  exists 
in  the  vicinity  of  a  planet,  it  is  necessary  to  consider  the  much  more 
complex  problem  of  "magnetic  field  plus  solid  body". 

It  was  Lemaitre  and  Vallarta  (1936a, b) ,  later  summarized  by  Vallarta 
(1961) ,  who  first  developed  a  solution  to  the  problem  in  terms  of  allowed, 
rather  than  forbidden,  access,  in  relation  to  the  geomagnetic  field.  Par  a 
given  rigidity,  they  realized  that  in  the  absence  of  the  Earth, 
trajectories  which  are  asymptotic  to  the  simplest  bound  periodic  orbits 
form  the  generators  of  a  cone,  within  which  all  possible  directions  of 
arrival  are  accessible  to  charged  particle  entering  the  field  from 
infinity.  Lemaitre  and  Vallarta  took  account  of  the  presence  of  the  earth 
by  introducing  as  additional  delimiters,  within  the  cane  of  allowed 
directions,  trajectories  which  are  at  seme  point  tangential  to  the  Earth's 
surface  and  whose  directions  of  approach  to  the  point  concerned  lie  within 
both  sets  of  delimiters  (Lemaitre  and  Vallarta,  1936b) ,  and  which, 
therefore,  are  allowed  under  both  considerations.  lemaitre  and  Vallarta 
also  used,  at  various  times,  the  terms  "main  cone"  and  "region  of  full 
light"  to  refer,  apparently  synonymously,  to  this  ccne. 

We  note  that,  unlike  the  Stormer  cone,  the  main  cone  was  initially 
defined  as  a  solid  quantity.  However,  as  in  the  case  of  the  Stormer  cone, 
but  rather  less  seriously,  seme  slightly  inconsistent  usages  have  appeared 
in  the  literature.  To  clarify  this  situation,  and  to  prepare  for  the 
symmetric  set  of  definitions  which  appears  in  section  5,  we  shall,  for  the 
remainder  of  this  report,  use  the  phrase  "allowed  cone"  to  describe  the 
solid  angle  region  of  full  light,  as  defined  by  Lemaitre  and  Vallarta,  and 
will  reserve  the  term  "Main  oone"  to  describe  the  boundary  of  the  alleged 
cone.  The  use  of  the  word  "oone"  can  be  misleading,  since  the  main  oone 
can  be  of  very  oonvoluted  shape,  and  can  even  break  into  two  or  more 
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separate  solid  angle  regions  (or  cones)  in  some  situations. 

The  allowed  ocne  is  a  very  useful  geometric  concept.  Par  any  given 
rigidity  and  location  in  a  magnetic  field,  a  solid  angle  region  can  be 
determined  within  which  cosmic  ray  particles  have  free  access  from  outside 
the  field.  For  decreasing  values  of  rigidity  the  size  of  the  allowed  cone 
diminishes,  until  at  seme  lower  value  it  vanishes. 

Lemaitre  and  Vallarta  (1936b)  used  a  mechanical  analogue  computer  to 
calculate  trajectories  by  a  numerical  integration  development  of  the 
initial  value  problem.  The  assumption  of  axial  symmetry  in  the  magnetic 
field  allowed  them  (like  Stormer  before  them)  to  restrict  their 
calculations  to  those  involved  with  tracing  particle  motion  in  the  meridian 
plane.  Use  of  the  Stormer  length  parameter 


(where  e  is  the  charge  of  the  particle,  H  is  the  dipole  magnetic  moment,  m 
is  the  particle  mass,  and  v  its  velocity)  permitted  generalization  of  their 
results  from  a  limited  number  of  traces  of  trajectories  asymptotic  to  bound 
periodic  orbits.  Having  thus  investigated  the  behavior  of  the  families  of 
asymptotic  trajectories,  Lemaitre  and  Vallarta  were  able  to  deduce  the  form 
of  the  allowed  cone  as  a  function  of  latitude  and  rigidity.  Later  work 
which  incorporated  quadripole  terms  did  not  alter  this  generality. 

When  the  question  of  cosmic  ray  access  to  a  specific  location  is 
considered  in  terms  of  the  effect  in  a  given  direction,  it  is  evident  that 
charged  particles  with  rigidities  in  excess  of  a  certain  limiting  value 
will  have  free  access  to  that  point.  This  limiting  rigidity  value  cor¬ 
responds  to  that  of  the  main  cone  in  the  direction  concerned  and  has  been 
termed  the  "main  cone  cutoff  rigidity4'  (Shea  et  al. ,  1965) .  A  particle 
arriving  in  this  direction  with  this  rigidity  value  will  have  traveled 
along  a  path  which  either  is  asymptotic  to  a  bound  periodic  orbit  or  which 
grazes  the  top  of  the  atmosphere  (in  the  case  of  the  Earth) .  The  Stormer 
cone  and  the  main  cone  can  never  overlap,  although  at  lew  latitudes  they 
can,  in  principle,  ocme  into  contact  over  a  limited  range  of  direction. 


Figure  1.1  shews,  conceptually,  the  spatial  relationship  between  the 
allowed  acne,  the  main  cone,  the  Stormer  cane,  the  forbidden  acne,  and  the 
cosmic  ray  penurrtora.  The  latter  is  the  solid  angle  zone  which  lies  between 
the  main  acne  and  the  Stormer  cane,  and  is  in  general  a  region  of 
alternating  allowed  and  forbidden  bands  of  directions  of  arrival. 

Cosmic  ray  access  to  the  penumbral  region  cannot  easily  be  evaluated, 
and  so  the  properties  of  the  terrestrial  penumbra  have  been  studied 
intensively  in  an  attempt  to  find  discemable  characteristics  by  which 
access  can  be  defined.  The  method  of  Lemaitre  and  Vallarta  was  applied  to 
the  study  of  the  larger  scale  forbidden  band  structures  in  the  terrestrial 
penumbra  by  later  investigators,  such  as  Hutner  (1939a, b) ;  Schremp 
(1938a, b) ;  Kasper  (1959)  and  Schwartz  (1959) .  Their  work  showed  the 
penumbra  to  possess  a  very  complex  structure  of  allowed  and  forbidden  bands 
of  varying  rigidity  widths.  Each  individual  forbidden  band  is  the 
manifestation  of  particular  lew  altitude  points  which,  for  a  given  family 
of  trajectories,  intersect  the  surface  of  the  Earth  (or  other  planet)  over 
a  finite  domain  of  rigidity.  It  was  not  until  the  advent  of  high  speed 
digital  computers  that  this  matter  could  be  further  investigated. 


1.4  CLASSICAL  TERMINOLOGY  AS  APPLIED  TO  DIGITAL  OCMEUIER  RESOLES 

The  digital  computer-based  studies  determine  the  access  of  cosmic  rays 
to  particular  locations  by  tracing  trajectories,  using  numerical 
integration,  through  the  planetary  field  as  represented  by  a  high  order 
mathematical  model.  As  the  arrival  location  is  given  at  the  cutset,  the 
trajectory  actually  traced  is  of  a  negatively  charged  particle  leaving  the 
location  in  the  direction  concerned.  This  trajectory  is  identical  to  that 
of  an  inbound  positive  particle  arriving  at  the  site  from  the  same 
direction.  If  a  particular  outbound  trajectory  is  unable  to  escape  from 
the  magnetic  field  (due  to  periodic  orbit  trapping,  or  due  to  the 
intersection  with  the  surface  of  the  planet)  then  clearly  a  particle  from 
infinity  is  unable  to  traverse  the  same  path  in  the  reverse  direction,  and 
this  arrival  direction  is  forbidden. 
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In  the  trajectory-tracing  technique,  the  form  of  the  allowed  and 
forbidden  access  regions  relating  to  any  point  in  the  field  may  be  sampled 
over  any  of  one  or  more  parameters  (rigidity*;  latitude,  longitude,  and 
altitude  of  location;  and  zenith  and  azimuth  angles  of  arrival)  by 
calculating  sets  of  trajectory  traces  with  incrementally  spaced  values  of 
the  given  parameter.  For  example,  where  it  is  desired  to  deduce  the 
periumbral  structure  in  a  given  direction  for  a  certain  location  as  a 
function  of  rigidity,  a  series  of  trajectories  would  be  calculated  at 
spaced  intervals  of  rigidity,  over  the  total  range  extending  from  a  value 
sufficiently  high  that  free  access  is  ensured,  through  to  a  low  value  at 
which  it  may  confidently  be  assumes  that  access  would  be  forbidden.  The 
size  of  the  rigidity  increment  chosen  depends  on  the  particular  requirement 
in  hand  -  relatively  large  increments  are  employed  where  only  gross  details 
of  the  penumbra  are  required,  and  fine  increments  if  more  detailed 
information  is  sought  (at  the  expense  of  a  proportionate  increase  in 
computer  run  time) . 

The  majority  of  recent  studies,  both  of  cut-offs  and  penumbral 
features  in  the  geomagnetic  field,  have  employed  rigidity  as  the  parameter, 
leading  automatically  to  considerations  of  particle  access  in  a  particular 
direction.  As  noted  in  section  1  we  intend  to  call  this  the  rigidity 
picture,  in  contrast  to  the  directional  picture  used  by  the  early 
investigators  in  their  studies  of  particle  access  as  a  function  of  arrival 
direction  ah  a  particular  rigidity.  This  and  the  following  section  discuss 
access  from  the  rigidity  viewpoint,  Whilst  section  5  provides  definitions 
appropriate  to  each  point  of  view. 

The  difference  between  the  reference  frames  used  in  the  earlier 
analytic  and  later  digital  cccputer-based  work  has  led  to  confusion.  Not 
having,  at  the  time,  a  full  realization  of  the  distinction  between  the  two 
viewpoints,  the  initial  users  of  the  digital  ccnputer  technique  applied  the 


*  Rigidity  -  momentum  per  unit  charge  is  a  more  convenient  unit  to  use  than 
is  energy,  since  the  results  obtained  are  then  independent  of  the  particle 
species. 
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1.5  PRESENT  UNDERSTANDING  OF  THE  FHENCMENOIOGY  OF  PARTICLE  ACCESS  TO 

POINTS  IN  A  PLANETARY  MAGNETIC  FIELD 

Recent  investigations  of  the  distributions  of  cutoff  values  at  various 
sites  and  altitudes,  and  examinations  of  the  penumbra  in  different 
situations,  as  well  as  development  of  more  sophisticated  techniques  for 
relating  particular  penumbral  features  to  specific  trajectory  forms,  have 
brought  a  better  understanding  of  the  phenomenology  of  charged  particle 
aocess  to  points  within  the  geomagnetic  field.  The  limited  ncn-terrestrial 
work  done  to  date  (the  examination  of  access  into  the  main  magnetic  field 
of  Jupiter,  Oooke  (1985) )  confirms  that  the  phenomenology  pertains  to  other 
planetary  fields,  as  expected.  In  the  light  of  this  current  understanding 
it  is  possible  to  examine  the  main  respects  in  which  the  earlier 
definitions,  developed  for  a  directional  view  of  access  in  a  field 
possessing  axial  symmetry,  fail  to  translate  to  the  rigidity  view  of  cosmic 
ray  access  to  points  in  the  real,  asymaetric,  field. 

Bound  periodic  orbits  have  especial  significance.  There  is  no  doubt 
that  they  can  exist  in  the  real,  asymmetric  field  -  the  existence  of  the 
terrestrial  trapped  radiation  belts  is  a  striking  confirmation  of  this 
fact.  Trajectory  calculations  suggest,  too,  that  quasi-bound  periodic 
orbits  can  be  expected  to  exist  in  the  real  field  of  both  Earth  and 
JUpiter.  It  is  necessary  to  investigate  the  significance  of  bound  periodic 
orbits  in  the  'real'  field,  to  ascertain  whether  or  not  orbits  of  the  kind 
associated  with  the  Stormer  cone  exist  in  this  field,  then  to  determine 
whether  access  to  the  point  concerned  from  outside  the  field  is  forbidden 
at  all  lower  rigidities.  It  is  normally  possible  to  infer,  by  calculation, 
that  a  direction  evidently  lies  within  the  cone,  by  observing  the  existence 
of  a  continuum  of  quasi-periodic  orbits  over  a  finite  range  of  initial 
parameters.  The  summary  plot  technique  developed  by  Oooke  and  Bredesan 
(1981)  in  relation  to  access  into  the  geomagnetic  field  represents  one  way 
of  detecting  such  a  continuum.  Because  of  the  departure  of  the  real  field 
from  axial  symmetry  it  could  not,  however,  be  assumed  that  entry  into  such 
orbits  is  unequivocally  forbidden.  (We  exclude  from  our  discussion  the 
time  variations  of  the  real  field  which,  of  course,  may  allow  entry  into 
otherwise  completely  bound  periodic  orbits  under  oertain  conditions. ) 
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In  regard  to  the  trajectory  forms  associated  with  the  main  cane, 
quasi-bound  periodic  orbits  have  particular  significance.  Although  this 
type  of  trajectory  is  a  singularity  (as  indeed  it  is  in  the  dipole  field, 
existing  only  for  a  single  value  of  energy) ,  it  can  be  identified,  because 
all  trajectories  within  a  finite  range  of  the  parameter  space  in  which  such 
an  orbit  is  located  show  a  characteristic  doubly  asymptotic  farm  (Iemaitxe 
and  Vallarta,  1936a) .  Figure  1.2  shews  a  real  field  near-main-oone 
trajectory  which  has  this  form.  It  is  possible  to  determine  a  value  of 
main  cone  cutoff  (or,  more  generally,  to  establish  the  position  of  the  main 
cone  in  any  parameter  spaoe)  by  ascertaining  the  trajectory  which  is 
characteristically  asymptotic  to  the  simplest  type  of  periodic  orbit.  We 
shall  refer  to  such  main  cutoffs,  which  are  purely  magnetically  controlled, 
as  being  of  type  I,  in  contrast  to  type  II  cutoffs  which  are  controlled  by 
trajectory-planet  intersections.  Because  the  radius  of  curvature  of  these 
orbits  is  large  relative  to  the  scale  dimensions  of  the  field,  the  question 
of  whether  the  orbits  are  truly  bound,  or  merely  possess  quasi-bound 
properties  over  large  trajectory  path  lengths,  has  proved  difficult  to 
resolve. 


In  order  to  illustrate  important  aspects  of  the  phenomenology  of 
charged  particle  approach  to  locations  within  a  planetary  magnetic  field, 
typical  characteristics  of  the  terrestrial  trajectory  shown  in  figure  1.3 
will  be  examined.  For  the  purpose  of  this  discussion  the  trajectory  will 
be  considered  as  possessing  two  portions,  identified  as  A  and  B  in  the 
figure.  In  portion  A,  a  positive  particle  enters  from  outside  the  field 
and  penetrates  downward  whilst  drifting  from  east  to  west  around  the  Earth 
towards  the  longitude  region  of  the  arrival  point.  It  is  in  this  part  of 
the  trajectory  that  hemisphere  to  hemisphere  oscillations,  associated  with 
the  quasi-bound  behavior,  are  to  be  found,  particularly  if  the  trajectory 
is  a  near-main-cone  one  (as  it  is  in  figure  1.3) .  The  characteristic 
quasi-bound  behavior  is  seen  in  a  more  exaggerated  form  in  the  trajectory 
of  figure  1.2.  Portion  B  of  the  trajectory  of  figure  1.3  describes  the 
final  motion  of  the  particle  towards  the  arrival  point.  Typically,  if  the 
arrived  point  is  at  other  than  very  low  latitudes  the  trajectory  will  loop 
one  or  more  times  within  the  local  field.  For  practical  purposes  the 
boundary  between  portions  A  and  B  of  the  trajectory  may  be  considered  to 
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FIGURE  1.2 

Typical  "real  field"  dcubly  asynptotic  trajectory  which  exists  at 
rigidity  values  close  to  the  main  cone  (trajectory  starting  point 
parameters:  geographic  latitude  10°  N,  longitude  270°  E,  zenith  angle 
7.9368°,  azimuth  270°,  rigidity  10.0  GV,  400  km  geocentric  altitude) 
calculated  using  the  International  Geomagnetic  Field  extrapolated  forward 
to  epoch  1980.0  (IAGA  Division  1  Study  Group,  1976).  This  trajectory  is 
referred  to  as  "doubly  asynptotic"  because  it  is  twice  asymptotic  to  the 
shell  of  essentially  constant  geomagnetic  radius  which  contains  the  central 
quasi-bound  section.  One  asyuptote  lies  in  the  trajectory  section  close  to 
its  entry  from  outside  the  field,  and  the  other  lies  close  to  the  final 
Earth  approach. 
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lie  at  the  extrema  latitudinal  excursion  before  the  equator  is  crossed  for 
the  last  time  as  the  particle  moves  towards  the  point  of  arrival.  There  is 
a  direct  correspondence  between  the  location  of  these  characteristic 
trajectory  forms  and  the  type  of  main  cone  produced.  In  particular,  type  I 
cutoffs  are  associated  with  quasi-bound  periodic  behavior  in  portion  A  of 
trajectories,  whilst  type  II  cutoffs  are  associated  with  earth  tangency  in 
trajectory  loops  which  lie  within  portion  B. 

In  general,  trajectory-planet  intersections  which  occur  in  portion  A 
produce  distinctly  different  structures  (penumbral  regions  inaccessible  to 
cosmic  rays)  than  do  those  which  occur  in  portion  B.  Consider  the  effect 
of  a  lew  point  intersecting  the  planet  within  portion  A.  Such  a  low  point 
can  only  occur  for  rigidities  at  or  below  the  main  cutoff,  because  only 
then  can  portion  A  possess  a  loop.  Above  the  main  cutoff  rigidity  value, 
the  trajectory  will  pass  once  through  the  "shell"  in  which  the  bound 
periodic  orbit  would  lie,  without  reversing  its  direction  of  motion 
relative  to  the  shell.  At  the  precise  main  cutoff  rigidity,  the  orbit  will 
lie  totally  within  the  shell.  At  rigidities  even  marginally  below  the  main 
cutoff,  the  trajectory  will  possess  reversals  in  the  sense  of  motion 
relative  to  the  shell,  resulting  in  the  formation  of  one  or  more  loops 
along  the  trajectory. 

The  rigidity  range  over  which  any  given  loop  will  exist  depends  upon 
the  path  distance  of  the  loop  from  the  final  point  of  arrival  in  the  field. 
If  ary  such  loop  intersects  the  surface  of  the  planet  then  a  penumbral 
forbidden  band  will  be  produced,  whose  width  in  rigidity  will  be  at  most 
equal  to  the  rigidity  range  over  which  the  intersecting  loop  exists.  At 
large  path  distances,  the  loops  and  therefore  ary  associated  penumbral 
bands,  can  only  exist  for  extremely  narrow  ranges  of  rigidity,  whereas  at 
short  path  distances  from  the  arrival  point  the  penumbral  bands  are 
relatively  stable  structures,  existing  over  rigidity  ranges  of  up  to  1  GV 
or  more. 


In  portion  B  of  the  trajectory,  particles  leave  the  quasi-bound 
section  of  the  orbit  and  commence  the  helical  motion  associated  with  their 
approach  to  the  final  arrival  point.  These  loops  exist  at  all  rigidities 
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and  (x2>i2)  i  then  a  and  7  are  given  by 
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Conversely,  if  the  location  of  one  point  is  known  (Ai,V>ir  say),  then 
it  is  possible  to  determine  the  latitude  and  longitude  of  a  second  whose 
position  is  defined  in  terms  of  a  specified  a  and  7 ,  by  means  of 


latitude  •  arcsin  (cos  0  sin  y  cos  Aj  +  cos  y  sin  Aj) 


cos  y  “  sin  A  sin  latitude 
longitude  -  ^  1  +  arcos  ( _ 1 _ 

cos  A  cos  latitude 
l 


III 


3)  The  calculation  of  the  magnetic  zenith  and  azimuth  uses  the 
strategy  of  setting  up  a  vector  of  known  length  fR")  pointing  in  the 
direction  of  interest.  The  position  of  the  tail  of  the  vector,  being  the 
location  of  interest,  is  of  course  known  in  both  geographic  and  magnetic 
coordinates  (the  latter  determined  by  step  1) ,  and  similarly  both 
geographic  and  offset  dipole  radius  values  pertaining  to  this  point  are 
known. 


The  geographic  coordinates  of  the  position  of  the  head  of  the  vector 
can  be  determined  by  using  calculated  a  and  7  values,  derived  as  follows 
(which  values  pertain  to  the  projection  of  the  vector  onto  a  spherical 
surface) . 


Y  -  srctsn  (R"  sin  te/(R  +•  R"  cos  re)) 
a  ■  at 


IV 


Having  evaluated  these  angles,  the  geographic  latitude  and  longitude  of  the 
vector  head  can  be  calculated  using  the  relationships  III.  The  offset 
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FIGURE  1.3 


The  trajectory  within  the  geomagnetic  field  of  a  5.65  GV  oomlc  ray 
particle  an  route  to  Williaastown,  Australia  (32.75°  S,  151.80°, 
geographic) ,  where  it  will  arrive  at  60°  zenith  angle  and  an  az  truth  of 
150°  east  of  geographic  north.  For  the  purposes  of  a  general  discussion  of 
the  effect  of  trajectory  characteristics  at  points  within  the  field,  two 
portions  of  this  typical  trajectory  are  identified,  and  are  labelled  as  A 
and 


up  to  values  significantly  above  the  type  I  main  cutoff,  and,  since  they 
lie  at  short  path  lengths  from  the  arrival  point,  are  stable  features  of 
the  trajectory.  However,  in  reality,  the  planet's  surface  does  prevent 
seme  of  the  trajectories  from  arriving  in  directions  which  might  otherwise 
be  allowed.  The  resulting  solid  angle  of  forbidden  directions,  at  any 
given  site,  is  referred  to  (by  analogy  with  optics)  as  the  planet's  cosmic 
ray  "shadow".  The  entire  range  of  directions  affected  (made  up  of  the 
shadows  associated  with  intersections  of  each  of  the  loops)  is  called  the 
"shadow  cone".  The  first  order  shadow  cane,  the  most  stable  of  these 
structures,  is  the  range  of  directions  associated  with  Earth  intersections 
of  the  loop  nearest  the  site;  the  seoend  order  shadow  cone  is  that 
associated  with  the  seoend  loop,  and  so  on.  When  the  shadow  cone  extends 
above  the  type  I  sain  cone  structure,  the  resulting  structure  is  easily 
recognizable.  Figure  1.4  illustrates  an  example  of  the  terrestrial  shadow 
oane  structure  as  shown  in  the  results  of  Cooke  and  Humble  (1970)  and 
Humble  and  Cooke  (1975) .  In  such  situations  it  is  the  shadow  cutoff 
structure  which  defines  the  sain  cone  over  the  range  of  directions 
affected. 

The  terrestrial  centered  dipole  field  shadow  cone  described  by  the 
generators  developed  by  Schremp  (1938a, b) ,  adopted  by  Vallarta  (1938)  and 
used  in  his  subsequent  work,  was  later  found  to  be  in  strong  disagreement 
with  the  experimental  measurements  of  Winckler  and  Andersen  (1954) ,  and  the 
later  trajectory  calculations  of  Schwartz  (1959)  and  Kasper  (1959) .  In 
addition,  it  should  be  noted  that  the  distinctions  between  shadow  and 
penumbral  forbidden  structures  discussed  by  Kasper  (1959,  I960) ,  Schwartz 
(1959) ,  and  Vallarta  (1961)  are  ocnpletely  equivalent  to  that  presented 
above.  The  earlier  interpretations,  however,  are  hard  to  follow  because  of 
the  difficulty  of  perceiving  all  dimensions  of  the  physical  reality 
associated  with  the  given  characteristics  of  trajectories  projected  onto 
the  meridian  plane.  Alpher  (1950)  converted  the  orthogonal  projections 
developed  by  the  earlier  workers  into  convenient  forms  to  enable  their 
comparison  with  experimental  data. 

In  addition  to  producing  the  shadow  cone,  trajectory  loops  within  the 
local  field  also  cause  systematic  structuring  of  the  directional 
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FIGURE  1.4 

Variation  of  upper  cutoff  with  geographic  azimuth  for  a  zenith  angle 
of  80 P  at  Brisbane,  Australia  (27.42°  S,  153.08°  E,  geographic).  The 
obtrusive  structure  between  80°  and  230°,  identified  by  the  hatched  lines, 
is  "Shadow  cone"  structure  associated  with  short  range  earth-intersecting 
trajectories  (after  Cooke  and  Rumble,  1970) . 
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distribution  of  type  I  cutoffs  at  nan  low-latitude  sites.  At  any  such  site 
the  main  cutoff  values  are  systematically  within  a  right 

circular  cane,  called  the  "locp  cone'1  by  Cooke  and  Humble  (1970) ,  the  axis 
of  which  lies  along  the  direction  of  the  local  field  at  the  site.  With 
increasing  latitude,  as  each  new  loop  forms  in  the  local  field,  a  new  cone 
develops,  and  increases  in  size. 

At  the  edge  of  each  loop  cone  the  main  cutoff  distribution  displays 
"folding’'  -  an  abrupt  change  in  cutoff  value.  Within  the  geomagnetic  field 
cutoff  changes  of  up  to  40%  have  been  found  in  association  with  the  folding 
at  the  edge  of  the  first  order  loop  cone,  the  acne  associated  with 
trajectories  having  a  single  loop  in  their  portion  B.  This  effect,  which 
is  visible  although  not  discussed  in  the  work  of  Lemaitre  and  Vallarta,  is 
particularly  obvious  in  azimuthal  scans  of  main  cutoffs  at  affected 
location.  Hie  folded  structure  can  be  such  that  a  curve  of  main  cutoffs 
plotted  for  a  range  of  azimuth  angles  at  constant  zenith,  can  be  seen  to 
turn  back  under  itself,  displaying  an  appreciable  overhang.  In  directions 
where  this  effect  occurs,  three  main  cones  exist  simultaneously,  one  far 
each  of  three  rigidities.  This  effect,  which  applies  only  over  a  limited 
range  of  directions,  produces  an  exception  to  the  otherwise  general 
assertion  that  in  any  direction  above  the  line-of -sight  horizon  the  main 
cutoff  defines  a  rigidity  limit  above  which  trajectories  of  all  rigidities 
are  allowed.  In  the  case  of  folded  structure,  only  the  highest  main  cutoff 
value  represents  such  a  rigidity  limit. 

Examination  of  the  form  of  individual  trajectories  reveals  general 
patterns,  such  as  path  length,  behavior  at  equator  crossings,  and  the 
development  of  characteristic  loop-like  structures,  which  relate  to 
distinct  fiducial  marks  recognizable  in  observations  of  the  cosmic  ray 
rigidity  spectrum  at  appropriate  locations.  For  example,  the  Initial 
development  of  loops  in  trajectories,  as  rigidity  is  reduced  at  the  sane 
site  and  arrival  direction,  signifies  that  the  rigidity  is  approaching  that 
of  the  sain  cone  in  that  direction.  Correspondingly,  the  "first 
discontinuity"  (section  5)  in  the  asyrrptotic  longitude  of  allowed 
trajectories  has  been  found  to  be  an  indicator  of  the  approach,  in  rigidity 
apace,  of  the  main  cutoff  (Shea  and  Smart,  1971;  Fludkiger  et  al. ,  1983). 


20 


the  then-existing  directional  picture  definitions  in  their  determinations 
of  cosmic  ray  cutoff  rigidities.  For  example,  the  highest  rigidity 
al lowed/ forbidden  trajectory  transition  found  in  any  given  direction  was 
labelled  as  the  main  cone  cutoff  (Shea  et  al. ,  1965) .  me  main  cone  must, 
of  necessity,  lie  at  a  rigidity  greater  than  or  equal  to  the  rigidity  value 
of  the  highest  detected  allowed/ forbidden  transition.  Ihe  highest 
forbidden  band  can  be  extremely  narrow  and  is  therefore  often  in 

actual  calculations  performed  ar  discrete  rigidity  increments  (Shea  and 
Smart,  1975) .  Thus  it  is  clear  that  the  highest  rigidity  al lcwed/ forbidden 
transition  confuted  in  a  systematic  survey  at  discrete  rigidity  intervals 
does  not  necessarily  identify  the  true  main  cone. 


Similarly,  use  of  directional  picture  definitions  in  a  rigidity  based 
investigation  led  Shea  et  al.  (1965)  to  identify  the  lowest  rigidity 
allowed/ forbidden  transition  found  as  the  Stormer  cone  cut-off,  and  so  this 
name  came  into  use  in  presenting  and  discussing  such  results  of  trajectory 
calculations,  even  though  it  has  since  become  clear  that  this  association 
in  most  situations  is  probably  invalid.  Vallarta  (1938)  stated  that  even 
in  a  staple  dipole  field  it  may  be  that  allowed  trajectories  do  not  lie 
adjacent  to  the  Stormer  cone.) 

Not  only  has  the  question  of  the  identification  of  the  "true"  main 
acne  and  Stormer  cutoffs  in  the  real  field  arisen  in  the  ocnputer-based 
studies,  some  additional  ocoplexities  have  also  become  apparent.  Far 
example,  investigations  by  Humble  et  al.  (1981)  of  the  pattern  of  access  to 
earth  satellites  have  shown  that  an  effectively  reversed  cutoff  situation 
exists  in  directions  below  the  optical  (line  of  sight)  horizon  far 
particles  approaching  with  relatively  high  rigidity.  Particles  of  infinite 
rigidity  are  unable  to  reach  a  satellite  from  below  horizontal  directions, 
since  their  path  is  blocked  by  the  presence  of  the  planet.  As  the  rigidity 
is  steadily  reduced,  a  value  is  eventually  reached  at  which  particles 
arriving  from  westerly  (for  the  earth)  directions  are  unable  to  reach  the 
satellite,  due  to  curvature  of  their  path  in  the  magnetic  field.  Humble  «t 
al.  (1981)  termed  this  cutoff  the  "allowed"  cutoff,  a  name  which  they  now 
believe  to  be  undesirable.  This  cutoff  is  defined  in  section  1.5  as  the 
"Horizon  Limited  Rigidity". 
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It  is  possible  to  utilize  such  characteristics  to  reduce  the  number  of 
trajectory  calculations  required  for  cutoff  determinations. 

The  form  of  the  penumbra  in  the  real  geomagnetic  field  has  also  been 
mapped  using  the  digital  ccmputer-based  trajectory  tracing  technique  (e.g. 
McCracken  et  al. ,  1962;  Freon  and  McCracken,  1962?  Shea  et  al. ,  1965;  Oooke 
and  Humble,  1979) .  Depending  of  the  nature  of  the  allowed/ forbidden 
structure  in  a  given  direction  the  penumbral  transmission  for  primary 
cosmic  rays  arriving  in  that  direction  at  locations  in  the  field  can  vary 
fran  zero  to  100%. 

Penumbral  trajectories  are  almost  invariably  considerably  more  complex 
than  are  those  lying  within  the  allowed  ccne,  and  mapping  by  use  of 
incrementally  spaced  parameters,  generally  rigidity,  consequently  requires 
large  amounts  of  computer  time.  In  a  manner  similar  to  that  indicated 
previously  in  connection  with  main  canes  it  is  possible  to  identify 
trajectory  characteristics  peculiar  to  many  penumbral  bands.  This 
systematic  behavior  of  penumbral  has  enabled  Cooke  (1982)  to  develop  the 
"trajectory  parameterization"  technique,  by  which  individual  penumbral 
structures  may  be  traced  by  relating,  in  an  automatic  ccrputer  procedure, 
the  structure  to  the  trajectory  feature  responsible  for  its  existence. 

(This  technique  also  forms  the  basis  of  an  automated  method  of  locating  the 
main  cutoff  to  high  precisian.) 

In  relation  to  access  into  the  earch's  field  the  most  stable  penumbral 
band,  that  which  is  associated  with  the  Earth  intersection  of  a  lev  point 
in  the  loop  which  lies  at  the  last  equatorial  crossing  in  the  quasi-bewnd 
periodic  portion  of  the  incoming  orbit  (normally  the  second-last  equatorial 
crossing  before  the  trajectory  reaches  the  arrival  point) ,  has  been  called 
the  "primary  band"  by  Lund  and  Sorgen  (1977) .  Because  of  the  offset  of  the 
Earth's  equivalent  dipole,  and  the  consequent  positioning  of  trajectory 
loops  relative  to  the  Earth's  surface,  the  primary  band  may  or  may  not, 
depending  upon  the  latitude  and  longitude  of  the  site,  be  visible  as  an 
isolated  forbidden  band  positioned  above  the  remainder  of  the  forbidden 
penumbral  band  structure. 


A  general  Indicator  of  the  stability  of  the  forbidden  penumbral  bands 
is  the  path  length  of  the  virtual  trajectory  (the  path  the  trajectory  would 
have  followed  had  the  planet  not  been  present)  fran  the  last  planet  inter¬ 
section  to  the  location  of  interest.  In  the  "reverse  direction  trajectory* 
tracing"  sense  this  corresponds  to  the  path  length  along  the  virtual 
trajectory  free  the  origin  of  calculation  to  the  first  intersection.  All 
relevant  terrestrial  cosmic  ray  studies,  free  both  the  directional  and  the 
rigidity  viewpoints,  support  the  attribution  of  large  scale  penumbral 
structure  (wide  penumbral  bands,  for  example)  to  trajectory  intersections 
with  the  Earth  at  relatively  short  path  distances  from  the  calculation 
origin,  and  that  of  fine  penumbral  structure  to  Earth  intersections  which 
occur  at  long  virtual  trajectory  path  lengths.  Although  the  forbidden 
structures  formed  by  the  short  range  virtual  trajectories  are  very  stable, 
changes  in  the  assumed  field  model  have  a  very  direct  influence  on  the 
penumbral  structure  formed  by  long  range  trajectory  path  lengths. 

Hie  question  arises  as  to  the  reality  of  fine  penumbral  structures  in 
a  planetary  magnetic  field.  In  the  case  of  the  earth  it  is  clear  that  the 
atmosphere  does  not  (contrary  to  the  assumption  made  by  most  investigators) 
act  as  a  simple  barriers  by  which  cosmic  rays,  having  penetrated  to  a 
certain  altitude,  are  cleanly  removed  from  the  incident  cosmic  ray  flux. 
Clearly  the  total  atmospheric  depth  traversed  by  the  charged  particle 
(which,  in  addition  to  altitude,  depends  directly  cn  the  trajectory 
curvature  and  configuration  in  the  grazing  section)  will  be  a  factor  in 
determining  whether  a  particle  will  survive  a  grazing  incidence  (see 
Lezniak  et.  al. ,  1975).  Also,  as  pointed  out  by  Petrou  and  Scutoul  (1977), 
the  atmospheric  depth  which  can  be  traversed  by  a  particle,  without  its 
loss,  is  a  function  of  the  identity  of  the  particle.  Thus  it  could  well  be 
that  fine  penumbral  terrestrial  structures  are  smeared  cut  and  essentially 
unobservable  experimentally. 

Time  variations  in  a  planetary  magnetic  field  caused  by  the 
fluctuating  pressure  of  the  solar  wind  and  by  the  rotation  of  the  planet 
with  respect  to  the  solar  wind  direction  are  an  additional  effect  which 
would  lead  to  the  smearing  cut  of  the  fine  structure. 
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In  brief  review,  it  can  be  stated  that  when  trajectories  calculated  in 
a  time-invariant,  but  otherwise  realistic,  model  of  the  Earth's  field  are 
examined,  all  the  analytically  identified  structures  identified  by  the 
earlier  workers  in  a  simple  axially  symmetric  field  can  be  recognized.  The 
limited  calculations  of  trajectories  in  a  higher  order  representation  of 
the  main  field  of  Jupiter  confirm  that  these  same  structures  are  found  in 
that  situation  too. 

There  is  good  reason  to  introduce  additional  terminology  in  order  to 
define  useful  quantities  naturally  associated  with  the  standard  sampling 
method  of  determining  the  real  field  cutoff  values.  These  terms  are 
presented  in  the  following  section,  together  with  the  "classical"  terms, 
suitably  qualified  to  allow  their  application  in  real  field  situations. 

1.6  DEFINITIONS 

The  definitions  referred  to  in  the  preceding  sections  are  listed  in 
table  1.1,  and  then  discussed  and  illustrated  later  in  the  section.  The 
terms,  as  listed,  are  subdivided  into  terms  referring  to  the  direction 
picture  and  the  rigidity  picture  (see  sections  1.5.1  and  1.5.2).  The  list 
is  not  exhaustive,  Exit  seeks  to  portray  the  most  useful  quantities  in  each 
situation. 

1.6.1  Directional  Definitions 

The  following  terms  are  appropriate  for  use  with  the  directional 
picture.  Each  definition  is  for  charged  particles  of  a  single  specified 
rigidity  arriving  at  a  particular  point  in  the  magnetic  field. 

ALLOWED  OONE:  The  solid  angle  containing  the  directions  of  arrival  of  all 
trajectories  which  do  not  intersect  the  planet  and  which  cannot  possess 
sections  asymptotic  to  bound  periodic  orbits  (because  the  rigidity  is  too 
high  to  permit  such  sections  to  exist  in  the  directions  of  arrival 
conocmed) . 


TABLE  1.1  Summary  of  terms  used  in  cutoff  calculations.  Quantities 
describing  phenomena  which  are  equivalent  in  the  two  pictures  are  listed  cn 
the  same  line. 


DIRECTION  PICTURE  RIGIDITY  PICTURE 


Allowed  Cone 
Main  Cone 

Shadow  cone 
Penumbra 
Penuataral  Band 

Sterner  Cone 
Forbidden  done 


Cutoff  Rigidity 

Main  cutoff  Rigidity 

First  Discontinuity  Rigidity 

Shadow  Cutoff  Rigidity 

Penumbra 

Penumbral  Band 

Primary  Band 

Stormer  Cutoff  Rigidity 

Upper  Cutoff  Rigidity 
lower  Cutoff  Rigidity 
Horizon  Limited  Rigidity 
Effective  Cutoff  Rigidity 
Estimated  Cutoff  Rigidity 


MAIN  CONE:  The  boundary  of  the  allowed  ocne.  The  main  cone  is  constituted 
in  part  by  trajectories  which  are  asynptotic  to  the  simplest  bound  periodic 
orbits  and  in  part  by  trajectories  which  graze  the  surface  of  the  planet. 
(In  the  case  of  the  Earth  the  "surface"  is  generally  taken  to  be  the  tqp  of 
the  effective  atmosphere.) 

PCRBUXEN  OCNE:  Ihe  solid  angle  region  within  which  all  directions  of 
arrival  correspond  to  trajectories  which,  in  the  absence  of  the  solid 
planet,  would  be  permanently  bound  in  the  magnetic  field.  Access  in  these 
directions  from  outside  the  field  is,  therefore,  impossible. 

STCSMER  OCNE:  The  boundary  of  the  forbidden  ocne.  m  an  axially  symnetric 
field  the  surface  forms  a  right  circular  ocne. 
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SHADOW  CCNE:  The  solid  angle  containing  an  directions  of  particle  arrival 
vhich  are  excluded  due  to  short  range  planetary  intersections  of  the 
approaching  trajectories  whilst  loops  within  the  local  field  line  bundle 
are  being  traversed. 

PENUMBRA:  The  solid  angle  region  contained  between  the  wain  ocne  and  the 
Stonier  ocne.  In  general  the  penumbra  contains  a  oonplex  structure  of 
allowed  and  forbidden  bands  of  arrival  directions. 

FENUMHRAL  BAND:  A  contiguous  set  of  directions  of  arrival,  within  the 
penumbra,  the  members  of  which  are  either  all  allowed  or  all  forbidden. 
Often,  within  ary  given  forbidden  bend  structure,  it  is  possible  to 
determine  that  a  number  of  individual  bunds,  each  attributable  to  the 
intersection  of  the  associated  trajectories  with  the  planet  in  a  different 
low  point,  are  overlapping  to  produce  the  entire  forbidden  band  structure 
observed. 


1.6.2  Rigidity  Picture  Definitions 

The  terms  defined  here  axe  used  in  relation  to  particles  arriving  at  a 
particular  site  within  the  planetary  magnetic  field  from  a  specified 
direction. 


COTOFF  RIGIDITY:  The  location  of  a  transition,  in  rigidity  space,  from 
allowed  to  forbidden  trajectories,  as  rigidity  is  decreased,  unless 
otherwise  defined  the  value  normally  quoted.  In  representing  the  results  of 
ocnputer  calculations  is,  for  practical  reasons,  the  rigidity  of  the 
allowed  member  of  the  appropriate  juxtaposed  allowed^  forbidden  pair  of 
trajectories  ocaputed  as  part  of  a  spaced  aeries  of  traces.  Sometimes  the 
term  is  employed  to  refer  to  the  location  of  a  notional  transition  from  one 
region  to  another,  for  exaaple,  at  the  Stormsr  ocna,  where  an  allOMd 
trajectory  may  not  perhaps  exist  at  all. 


MAIN  COTOFF  RZGIDZTY,  %:  The  rigidity  value  at  which  the  direction 
concerned  is  a  generator  of  the  main  oone  as  defined  in  the  "direction 


picture".  The  associated  trajectory  is  either  one  which  is  asymptotic  to 
the  simplest  bound  periodic  orbit,  or  (owing  to  the  presence  of  the  solid 
planet}  is  one  which  is  tangential  to  the  planet's  surface. 


FIRST  DISQCNTINUITY  RIGIDITY,  R^  The  rigidity  associated  with  the  first 
discontinuity  in  asymptotic  longitude  as  the  trajectory  calculations  axe 
performed  for  successively  lower  rigidities,  starting  with  sane  value 
within  the  allowed  cone.  The  value  of  R]_  is  approximately  equal  to  the 
main  cone  cutoff  as  defined  above. 

SHADOW  CUTOFF  RIGIDITY,  The  rigidity  value  at  which  the  edge  of  the 

shadow  cane  lies  in  the  direction  concerned. 

PENUMBRA:  The  rigidity  range  lying  between  the  main  and  the  Stormer  cutoff 
rigidities. 


FENUMHRAL  BAND:  A  continuous  set  of  rigidity  values,  within  the  penumbra, 
all  members  of  which  have  the  same  general  access  characteristics,  either 
all  allowed  or  all  forbidden.  Often,  within  any  given  forbidden  band 
structure,  it  is  possible  to  determine  that  a  number  of  individual  bands, 
each  attributable  to  the  intersection  with  the  solid  planet  of  the 
associated  trajectories  in  a  different  low  point,  are  overlapping  to 
produce  the  entire  forbidden  band  structure  observed. 

PRIMARY  BAND:  The  stable  forbidden  penumbral  band  which  as  associated  with 
the  planetary  intersection  of  a  low  point  in  the  loop  which  lies  at  the 
last  equatorial  crossing  before  the  trajectory  (of  its  virtual  extension  in 
the  assumed  absence  of  the  planet)  takes  on  the  characteristic  guiding 
center  motion  down  the  local  field  line  bundle. 

STORMER  CUTOFF  RIGIDITY,  Rg:  The  rigidity  value  for  which  the  Stormer  cone 
lies  in  the  given  direction.  In  a  dipole  field  (and  perhaps  also  in 
real  geomagnetic  and  other  planetary  fields)  aooess  for  particles  of 
rigidity  values  lower  than  the  Stormer  cutoff  rigidity  is  forbidden 
outside  the  field.  In  a  dipole  approximation  to  the  magnetic  field, 
form  of  the  Stormer  equation  gives  the  Stormer  cutoff  rigidity,  in  GV,  as: 
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Rs 

r2[l  +  (1  -  cos3*.  cosE  $1nZ)1/2j2 

where  M  is  the  dipole  moment,  which,  far  the  geomagnetic  field,  has  a 
normalized  value  of  59.6  when  r  is  expressed  in  units  of  earth  radii,  x  is 
the  magnetic  latitude,  r  is  the  distance  from  the  dipole  in  earth  radii,  E 
is  the  azimuthal  angle  measured  clockwise  fran  the  geomagnetic  east 
direction  (far  positive  particles) ,  and  Z  is  the  angle  from  the  local 
magnetic  zenith  direction. 

UPPER  cutoff  RIGIDITY,  %:  The  rigidity  value  of  the  highest  detected 
allowed/forbidden  transition  among  a  set  of  ccnputed  trajectories.  The 
upper  cutoff  rigidity  can  correspond  to  the  main  cutoff  if  and  only  if  no 
trajectories  asynptotic  to  bound  periodic  orbits  lie  at  rigidities  higher 
than  this  value.  This  can  be  identified  from  the  nature  of  the  trajectory 
associated  with  the  main  cutoff. 

LOWER  cutoff  RIGIDITY,  The  lowest  detected  cutoff  value  (i.e.  the 
rigidity  value  of  the  lowest  allc*«y  forbidden  transition  observed  in  a  set 
of  ocnputer  calculations) .  If  no  penuntoa  exists,  will  equal  Ru- 

HORIZON  LIMITED  RIGIDITY,  The  rigidity  value  of  the  most  rigid  allowed 
trajectory  found  in  a  set  of  computer  calculations  performed  far  e  below 
horizon  direction  at  a  location  above  the  surface  of  the  planet. 

EFFECTIVE  CUTOFF  RIGIDITY,  R c?  The  total  effect  of  the  penuobral  structure 
in  a  given  direction  may  usefully  be  represented,  for  many  purposes,  fey  the 
“effective  cutoff  rigidity”  -  a  single  numeric  value  which  specifies  the 
equivalent  total  accessible  cosmic  radiation  within  the  penumbra  in  a 
specific  direction.  “Effective  cutoffs"  may  either  be  linear  averages  of 
the  allowed  rigidity  intervals  in  the  panumbra  (Shea  et  si.,  1965),  or 
functions  weighted  for  the  cosmic  ray  spectrum  and/or  detector  response 
(Shea  and  smart,  1970;  Dorman  et  el.,  1972).  Far  e  linear  weighting  this 
would  have  the  form: 


Rc  ■  R“  -  l  4R<  (allowed)  ; 

where  the  trajectory  calculations  were  performed  at  rigidity  intervals  R^. 

ESTIMATED  CUTOFF  RIGIDITY,  Rest:  A  value  obtained  by  using  an  empirically 
normalized  equation  to  approximate  the  cosmic  ray  cutoff  variation  in  the 
location  of  a  particular  point  within  a  magnetic  field  in  order  to  estimate 
a  cutoff  pertaining  to  the  point.  This  value  can  be  found  by  use  of  a 
variety  of  interpolation  techniques,  one  of  which  is  application  of  the 
Starmer  equation  given  previously.  Because  the  Stormer  equation 
characterizes  the  spatial  variation  of  the  cutoff  rigidity,  with 
appropriate  normalization  it  may  be  used  to  obtain  useful  estimates  of  the 
various  cosmic  ray  cutoff  rigidities,  over  intervals  of  latitude, 
longitude,  zenith  and  azimuth,  for  example.  In  practice,  estimates  of  the 
value  of  any  cutoff  can  be  obtained  from  adjacent  calculated  values  to  a 
reasonable  accuracy  by  employing  this  method. 


1.7  niUSTRATIGN  OF  DEFTNITTCNS 

Not  all  the  cutoffs  defined  in  the  previous  section  exist  far  every 
location  and  direction.  Figure  1.5  shows  typical  cutoffs  obtained,  using 
trajectory  calculations,  for  an  arbitrarily  selected  location  on  the  Earth 
(20°  N,  270°  E,  altitude  400  km)  at  a  zenith  angle  of  60°.  Forbidden 
rigidities  are  shown  in  black,  and  allowed  rigidities  are  indicated  by 
white.  The  upper  cutoff  rigidity  (Rq)  and  the  lower  cutoff  rigidity  (Rl) 
far  each  azimuth  can  readily  be  identified.  In  sene  directions,  azimuth 
angles  of  105°  and  255°  in  this  instance,  %  and  Rl  have  the  same  value, 
within  the  1%  rigidity  intervals  used  for  these  calculations.  Note  that  Ru 
and  Rl  will  be  equal  in  any  direction  far  which  no  allowed  penunfaral  bands 
exist,  a  situation  which  is  particularly  oawnon  at  lew  latitude  locations. 

Typical  cosmic  ray  aooess  patterns  at  zenith  angles  larger  than  that 
of  the  Earth  horizon  seen  fron  a  satellite  are  illustrated  in  figure  1.6, 


29 


RIGIDITY(GV) 


LATITUDE  N  20°  ALTITUDE  400km 

LONGITUDE  E  270°  ZENITH  60° 

,2.o  , - 


AZIMUTH  ANGLE  (CLOCKWISE  FROM  NORTH) 


FIGURE  1.5 

Illustration  of  cutoffs  obtained  at  a  particular  location  (20°  N,  270° 
E,  400  Job  altitude)  at  a  zenith  angle  of  60°.  Trajectories  at  1%  rigidity 
•pacing  intervals  were  calculated  every  30°  in  azlsuth  angle.  Ttorbldtat 
trajectories  are  illustrated  by  the  dark  areas,  and  allowed  trajectories  by 
white. 
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for  the  same  location  as  in  figure  1.5.  The  zenith  angle  of  120°  is  below 
the  direction  to  the  Earth's  horizon  (approximately  109°  far  a  spacecraft 
at  an  altitude  of  400  km) .  For  this  location  and  zenith  angle  all 
rigidities  are  forbidden  in  the  north  and  easterly  directions.  However, 
for  southerly  and  westerly  directions  a  range  of  allowed  rigidities  exists 
within  which  cosmic  ray  particles  may  approach  free  angles  below  the 
satellite/Earth  horizon  without  encountering  the  Earth.  In  such 
directions,  of  course,  extremely  rigid  cosmic  ray  particle  trajectories 
will  intersect  the  Earth  and  the  particles  are  therefore  unable  to  reach 
the  site.  The  highest  rigidity  that  has  access  at  any  specified  zenith 
angle  belcw  the  spacecraft/Earth  horizon  is  referred  to  as  the  horizon 
limited  rigidity.  These  horizon  limited  rigidities  are  evident  at 
azimuthal  angles  of  165°  through  315°  in  figure  1.6. 

Figure  1.7  illustrates  penumbral  structure  encountered  in  the 
determination  of  the  effective  vertical  cutoff.  He,  for  three  North 
American  locations  using  linear  weighting  of  penunhral  bands.  Note  that 
there  is  an  inherent  assumption  involved  in  determining  the  values  of  Be, 
that  the  result  of  a  trajectory  calculation  at  a  very  specific  rigidity 
typifies  the  result  over  a  discrete  rigidity  interval.  The  allowed 
intervals  in  the  penumbra  are  then  summed  (with  appropriate  weighting  far 
detector  response  and  spectrum  as  required)  and  subtracted  from  the 
calculated  upper  cutoff  rigidity. 


1.8  CONCLUDING  REMARKS 

It  should  be  borne  in  mind  that,  because  the  definitions  have 
deliberately  been  kept  usefully  general,  the  application  of  the  terms  may 
require  more  detailed  qualification  in  seme  individual  circumstances.  For 
example,  in  below  line-of -sight  horizon  directions  at  satellite  altitudes 
it  may  be  that  there  is  not  always  aocess  into  the  allowed  oona  for 
particles  of  any  given  rigidity,  and  so  consequently  there  may  not  be  a 
main  cutoff  in  seme  cases.  Nevertheless,  particles  within  a  limited 
rigidity  range  often  can  approach  and  have  aocess  in  such  below  hjrizon 
directions  (figure  1.5),  and  characteristic  aHmecyforbidden  banded 
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AUTITUOE  400  km 


LONGITUDE  £  270* 


o  80 
a 

TO 


45  75  05  135  (65  195  22S  255  295  515  345 

azimuth  angle  (Clockwise  from  north) 


FIGX23E  !•( 


Illustration  of  "allowed"  rigiditiaa  at  a  zenith  angle  balov  ths 
direction  of  the  local  earth  horizon  far  a  satellite  at  400  Xe  altitude  at 
a  location  20 P  H,  270°  E.  ftabiddan  trajectories  are  indicated  by  black 
and  allowed  trajectories  by  white.  The  "horizon  Halted  rigiditiaa"  are 
apparent  as  the  high  rigidity  transition  free  allowed  to  farbidtai 
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MT.  WASHINGTON  {44.30N,  288.70E)  R--I.43GV 


i.ro  i~20  rto  rio  Tso  1.6O  '  ’ 

RIGIDITY  (GV) 


NEWARK  (39.68N,  284.25E)  RC«2.I0  GV 
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RIGIDITY  (GV) 


CLIMAX  (39.37N,  253. 82 E)  RC-2.95GV 


FIGURE  1.7 


Illustration  of  the  relationship  between  the  penunferal  structure  and 
effective  vertical  cutoff  rigidity,  Rq,  for  three  North  American  cosmic 
«tati«w.  For  these  calculations  the  trajectories  were  vertically 
idenfc  to  the  geoid  at  the  station  location,  and  were  at 

rigidity  intervals  of  0.01  GV  using  the  International  Geomagnetic  Reference 
Field  for  1980.0  (Peddle,  1982) .  Forbidden  direction  of  arrival  are 
illustrated  by  the  dark  areas  and  allowed  by  white. 
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penumbral  structures  clearly  can  exist  in  these  situations.  The  question 
then  arises  as  to  what  is  the  meaning  of  the  term  penumbra  in  these 
directions.  With  suitable  qualification  (i.e.  of  the  rigidity  or  other 
range  considered)  the  ♦’penumbra"  can  be  defined  far  any  particular 
situation.  For  example,  the  "calculated  penumbra"  could  be  defined  as  the 
rigidity  interval  between  the  calculated  upper  cutoff  rigidity  and  the 
calculated  lowest  cutoff  rigidity. 

Experimental  cutoffs  are  beccndng  sufficiently  precise  to  allow  the 
resolution  of  the  gross  characteristics  of  the  terrestrial  cosmic  ray 
pertunfetra  (Lund  and  Sorgen,  1977;  Byrnak  et  al. ,  1981;  and  Soutoul  et  al. , 
1981) .  For  cutoff  values  derived  from  experimental  measurements, 
definitions  comparable  to  those  discussed  in  the  previous  section  can  be 
used.  In  this  situation,  the  terms  "upper  cutoff  rigidity"  and  "lower 
cutoff  rigidity"  should  be  prefaced  by  an  appropriate  qualifier  such  as 
"measured"  or  "experimental". 

There  is  no  doubt  that  other  physically  meaningful  quantities  exist. 

It  is  believed,  however,  that  the  cutoff  concepts  described  in  this  chapter 
presently  have  the  greatest  significance,  and  that  the  use  of  these 
definitions  should  alleviate  most  of  the  existing  confusion,  and  satisfy 
the  current  requirements  of  investigators  involved  in  cosmic  ray  amass 
studies,  both  in  relation  to  the  Earth  and  the  other  planets  possessing 
significant  magnetospheres. 


CHAPTER  2:  USE  OF  THE  STCSMER  EQUATION  TO  REPRESENT  DIRECTIONAL  CUTOFFS 


2.1  INTRODUCTION 

Far  many  purposes  It  is  useful  to  have  a  fast  means  of  estimating 
directional  cosmic  ray  cutoffs  pertaining  to  ary  specified  direction  and 
location.  The  estimation  of  "real"  directional  cutoffs  by  ocnputer  is  a 
lengthy  process,  expensive  in  ocnputer  tine  (even  if  the  techniques 
described  later  in  this  report  dramatically  increase  the  efficiency  of  this 
process) .  The  Sterner  cutoff  function  (Stormer,  1930,  1955) ,  which 
expresses  the  dependence  of  the  Stormer  cutoff  cn  location  and  direction  in 
a  dipole  approximation  to  the  Earth's  field,  offers  a  means  of  determining 
cutoffs  which  is  many  orders  of  magnitude  faster  than  ary  trajectory 
tracing  method,  but  one  which  is  sufficiently  precise  to  produce  useful 
estimates  in  non-critical  applications  -  far  the  Earth's  field  and  far  any 
other  planetary  field  which  may  be  approximated  by  a  dipole. 

The  fundamental  inprecision  in  the  Stormer  expression  in  representing 
real  field  cutoffs,  in  particular  due  to  failure  to  take  into  account 
higher  order  field  harmonic  terms,  or  the  width  and  transparency  of  the 
penumbra,  is  normally  exacerbated  by  the  use  of  centered  coordinates  when 
invoking  the  expression.  It  is  possible  to  appreciably  lsprove  the 
accuracy  of  the  cutoff  estimates  by  using  "magnetic"  coordinates  (l.e. 
offset  dipole  latitude,  longitude,  zenith,  and  azimuth)  when  employing  the 
expression,  and  in  this  way  to  take  into  account  the  offset  and  tilt  of  the 
equivalent  dipole.  By  this  means,  inherently,  the  effect  of  ignoring  the 
higher  order  field  terms  is  minimized.  Smart  and  Shea,  1977,  have 
discussed  the  advantages  of  using  offset  dipole  coordinates  in  conjunction 
with  the  Stormer  expression,  and  the  use  of  this  expression  for 
interpolating  cosmic  ray  cutoffs  over  intervals  within  which  precise 
calculated  values  do  not  exist. 

A  self  contained  system  for  transforming  from  geographic  coordinates 
to  offset  dipole  coordinates  is  described  here.  It  has  been  used  in  the 
calculation  of  cutoff  distribution  functions  for  use  in  estimating  detector 
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response  to  neutrinos  of  atmospheric  origin.  It  has  been  necessary  to  take 
into  account  the  effects  of  the  magnetic  cutoffs  over  the  entire  solid 
angle  of  acceptance  of  detectors  at  locations  within  the  geomagnetic  field. 


2.1  STCSMER  EQUATION  IMPLEMENTATION 

The  coordinate  transformation  has  a  number  of  etages,  which  are 
individually  described  in  the  following  (this  process  is  described  as  it 
applies  to  the  Earth's  field) : 

1)  A  coordinate  conversion  is  used  to  determine  the  offset  dipole 
latitude,  longitude,  and  radius  tram  the  nominated  geographic  latitude  and 
longitude,  and  geocentric  altitude.  This  conversion  takes  into  account  the 
offset  and  inclination  of  the  Earth's  equivalent  dipole  far  any  required 
epoch,  and  assumes  that  the  Earth  is  an  oblate  spheroid  of  eccentricity 
0.00674.  The  angle  conversion  equations  are  as  follows: 


offset  dipole  longitude  +  •  arctan  (R  sin  f  cos  1  -  y  cos 

C 

offset  dipole  latitude  9  *  arctan  (cos  4  tan  a)  I 


C 

offset  dipole  radius  R'  - 

cos  6  cos  4 


where  a  ■  arctan  (F/G)  +  a 


C  ■  C  cos  a 
cot(a  -  a) 


and  F  •  R  sin  X  -  x 

G  •  y  ainb  +  R  cot  X  cos  ( 9  -  c) 


x  and  i>  are  the  geographic  latitude  and  longitude  respectively;  R  is  the 
geographic  radius  at  the  specified  location;  x  is  the  displacement  of  the 
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dipole  from  the  center  of  the  earth  in  a  direction  parallel  with  the 
geographic  equatorial  plane;  a  is  the  inclination  of  the  dipole  axis  from  a 
direction  parallel  to  the  geographic  N-S  axis;  b  is  related  to  the  angle 
between  the  zero  magnetic  longitude  and  the  geographic  longitudinal  towards 
which  the  dipole  is  displaced?  and  c  is  the  geographic  longitudinal 
direction  parallel  to  which  the  direction  zero  offset  dipole  longitude 
lies. 


A  procedure  for  establishing  the  position  of  the  equivalent  dipole 
with  respect  to  the  geocentric  coordinate  system  at  any  epoch  is  presented 
by  Smart  and  Shea,  1977.  This  procedure,  which  makes  use  of  the  low  order 
terms  in  the  spherical  harmonic  representation  of  the  geomagnetic  field,  is 
discussed  more  fully  by  Roederer,  1972;  and  Akasofu  and  Chapman,  1975.  The 
values  of  x,  y,  a,  b,  and  c  used  in  the  presently  reported  coordinate 
transformation  can  thus  be  determined  as  required.  The  equivalent  dipole 
position  determination  is  included  as  a  subroutine  in  the  ocnputer 
which  executes  the  Stanner  cutoff  evaluation. 

2)  The  following  steps  rely  on  the  use  of  a  particular  means  at 
specifying  the  relative  position  of  two  points  on  the  surface  of  a  sphere. 
In  particular,  two  angles  are  used,  one  (a)  defines  the  angle  between  the 
great  circle  connecting  the  two  points  and  the  meridian  line  intersecting 
one  of  them,  and  the  other  (7)  is  the  angle  subtended  at  the  center  of  the 
sphere  by  the  two  points  (see  figure  2.1) . 

Figure  2.1  Diagram  defining  the  angles 
a  and  7  used  to  express  the 
relative  position  of  the  two 
points  Pi  and  Pj  on  a  spheri¬ 
cal  surface.  C  is  the  oenter 
of  the  sphere. 


If  the  two  points  are,  as  specified  by  latitude  and  longitude,  (>1,71) 
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dipole  1:  uda  and  longitude,  and  the  distance  R*  of  the  vector  head  from 
the  dipole  are  then  calculated  using  relationships  I.  Now,  having  the 
offset  dipole  coordinates  of  both  ends  of  the  vector,  and  angles  and 
relating  the  two  points  (in  the  offset  dipole  frame  of  reference)  can  be 
calculated  by  means  of  the  relationships  H.  The  azimuth  of  the  vector  in 
the  offset  dipole  frame  of  reference  is  simply  the  a  value,  whilst  the 
zenith  angle  is  given  by 

ze  “  ercoz  (R*  cos  *  -  R* ) 


where  R*  is  the  offset  dipole  radius  of  the  vector  tail  (i,e,  the  distance 
of  the  vector  tail  from  the  offset  dipole  center) . 


2.3  CALCUIATICN  OF  STCRMER  CUTOFF  VAIXJE: 


By  using  the  coordinate  transformation  procedure  described,  the 
position  (latitude,  longitude,  and  radius)  of  the  site  location  relative  to 
the  offset  dipole,  and  the  zenith  and  azimuth  pertaining  to  the  direction 
of  interest,  may  be  calculated  from  the  specified  geographic  coordinates 
and  geocentric  altitude.  At  each  step  during  the  derivation  of  these 
parameters  tests  are  performed  to  ensure  that  calculated  angle  values  lie 
in  the  correct  quadrant,  and  appropriate  corrections  are  made  if  they  do 
not.  Having  thus  determined  the  angles  and  distance  relative  to  the 
equivalent  dipole  centered  frame  of  reference  the  Starmar  equation  can  be 
invoked  with  the  greatest  possible  precisian. 


The  appropriately  normalized  Starmar  expression  is  as  follows: 


cutoff  • 


4 

_ 59.4  coz  9 _ 

r*2  (i  +  / 1  -  coz^  9  sin  as  ain  a#)* 


This  expression  takes  in  offset  dipole  latitude  (f),  zenith  (ze),  azimuth 
(az  -  measured  clockwise  from  magnetic  north) ,  and  redial  distance  from  the 
effective  dipole  center  (R*) ;  and  produces  cutoff  rigidity  values,  in  units 
of  GV. 
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2.4  USE  OF  STCSMER  CUTOFFS  IN  THE  EERIVATTCN  OF  GEOMAGNETIC  CUTOFF 

DISTKEBUTTCN  FUNdTCNS  PCR  USE  IN  ESTIMATING  EETECICR  RESPONSE  TO 

NEUHUNCS  OF  ATMOSPHERIC  ORIGIN 

A  procedure  has  been  developed,  utilizing  the  Sterner  cutoff 
expression  far  deriving  functions  which  characterize  the  effect  of 
geomagnetic  cutoffs  on  the  charged  primary  cosmic  rays  that  give  rise  to 
neutrinos  arriving  in  any  given  direction  at  specified  points  on  or  in  the 
earth.  These  cutoff  distribution  functions,  for  vise  in  atmospheric 
neutrino  flux  calculations,  have  been  determined  for  eight  nucleon  decay 
experiment  sites,  vising  a  technique  which,  in  addition  to  enplcying  the 
Stormer  expression,  assumes  col  inear  motion  of  neutrino  and  parent  primary. 

Some  large  cosmic  ray  detectors,  such  as  those  used  in  nucleon  decay 
research,  have  a  finite  response  to  neutrinos.  It  is  necessary  to  be  able 
to  calculate  the  expected  background  flux  of  neutrinos  incident  upon  such 
detectors  on  order  to  establish  whether  the  angular  distribution  of  certain 
classes  of  events  has  the  characteristics  of  nucleon  decay  or  of  neutrino 
interactions,  a  distinction  which  has  important  physical  and  astrcphysical 
implications. 

The  neutrinos  produced  in  the  atmosphere  (muon  and  electron  neutrinos 
associated  with  xtuon,  picn  and  kaan  decay)  constitute  a  major  part  of  the 
background.  The  rate  of  production  of  these  neutrinos  is  related  to  the 
intensity  of  primary  cosmic  rays  incident  upon  the  atmosphere,  which  in 
turn  depends  on  the  primary  cosmic  ray  spectrum,  and,  as  discussed  by 
Gaisser  (1982) ,  on  the  geomagnetic  cutoffs  pertaining  in  the  given 
situation. 


Let  the  differential  intensity  of  neutrinos  reaching  a  detector  in  the 
presence  of  the  geomagnetic  field  be  A(R,5 ,^) ,  where  R  is  rigidity,  and  9 
and  are  the  zenith  and  azimuth  angles  of  neutrino  arrival  at  the 
detector.  If  B(R,*,*)  is  the  intensity  of  neutrinos  that  would  exist  in 
the  absence  of  the  field,  then  the  functions  A  and  B  can  be  related  as 
follows: 

A(R,*,*)  •  B(R,#,y>)  C(R,#,*)  2.1 
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where  C(R,$  ,$)  expresses  the  effect  of  the  geomagnetic  cutoffs.  This 
"cutoff  distribution  functicn"  essentially  describes,  far  a  detector  in  a 
given  location,  the  fraction  of  the  total  solid  angle  of  the  detector 
accessible  to  neutrinos  that  are  descended  fran  charged  primaries  of  any 
given  rigidity. 

In  deriving  these  functions,  the  cutoffs  relating  to  points  over  the 
entire  earth's  surface  have  to  be  taken  into  account.  (Because  neutrinos 
can  penetrate  the  earth,  their  production  in  the  atmosphere  on  the  distant 
side  of  the  Earth,  as  veil  as  in  the  local  atmospheric  mass,  is 
significant.) 

A  ocnputer  based  procedure  has  been  developed  for  calculating  the 
functions.  It  assumes  that  the  detector  has  an  isotropic  directional 
response  to  neutrinos.  The  entire  4*  steradian  field  of  view  is  divided 
into  a  number  of  equal  solid  angle  zones,  where  each  zone,  essentially  an 
annulus,  lies  between  defined  upper  and  lower  zenith  angle  limits.  A 
separate  function  is  determined  far  each  of  the  zones.  (In  the  present 
analysis  eight  */2  steradian  zones  axe  employed,  whose  zenith  angle  limits 
are  listed  in  Table  2.1.) 

As  a  first  step  towards  deriving  the  functions,  each  solid  angle  zone 
is  divided  into  a  large  number  of  smaller  elements  of  equal  solid  angle. 

The  axial  zenith  and  azimuth  angles  are  confuted  far  each  element,  and  then 
these  angles,  together  with  the  specified  site  position,  axs  processed  to 
determine  the  location  and  angle  at  which  the  "  1  ine-of-sight"  vector 
intersects  the  assumed  production  level  in  the  atmosphere  (a  height  of  20 
km  above  the  earth,  in  the  reported  calculations) . 

A  further  set  of  transformations  converts  these  angles  into  magnetic 
coordinates.  In  this  work  and  inclined,  offset,  magnetic  dipole 
approximating  the  1980  magnetic  field  is  assumed.  The  angle  conversions 
and  cutoff  expressions  from  the  preceding  section  are  used  to  produce  the 
directional  cutoff  values. 
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TABLE  2.1 


Zenith  angle  extent  of  the  eight  solid  angle  zones  far  which  cutoff 
distribution  functions  have  been  calculated  at  each  site.  (The  zenith 
angle  is  expressed  relative  to  the  detector.  Thus  the  first  rone, 
extending  from  0°  to  41.4°,  accepts  neutrinos  traveling  downward,  whilst 
zone  8,  extending  from  138.6°  to  180.0°,  accepts  neutrinos  traveling 
upwards.) 

Zone  Zenith  Angle  Range  Weighted  Mean  Zenith  Angle 


(degrees) 

(degrees) 

1 

0.0 

- 

41.4 

27.7 

2 

41.4 

- 

60.0 

51.1 

3 

60.0 

- 

75.0 

67.9 

4 

75.5 

- 

90.0 

82.8 

5 

90.0 

- 

104.5 

97.2 

6 

104.5 

- 

120.0 

112.1 

7 

120.0 

- 

138.6 

128.9 

8 

138.6 

- 

180.0 

152.3 

These  cutoff  values  do  not  take  into  account  the  ncn-dipole 
components  of  the  Earth's  field,  penunhral  effects,  or  the  Earth's  cosmic 
ray  "shadow".  Far  these  reasons  there  are  significant  disparities  between 
calculated  real  field  cutoffs  and  the  Stormar  estimates.  It  is  estimated, 
on  the  basis  of  direct  comparative  calculations,  as  well  as  by  checks  mads 
with  the  use  of  the  published  real  field  cutoffs  of  Shea  and  smart  (1982) , 
that  the  agreement  is  within  about  40%  in  directions  where  the  shadow 
effect  is  not  present.  At  zenith  angles  greater  than  about  60°,  in 
directions  where  the  shadow  effect  in  present,  discrepancies  of  the  order 
of  a  factor  of  two  can  be  encountered.  Nevertheless  it  is  felt  that.  In 
this  first  approach,  the  S termer  cutoff  values  constitute  acceptable 
approximations  to  the  real  field  cutoffs. 

Having  thus  the  means  for  quickly  estimating  the  required  directional 
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cutoffs,  the  summing  of  the  effects  of  the  cutoffs  over  all  the  solid  angle 
elements  within  a  zone  is  carried  out  by  numerical  integration,  as  follows: 


C(K.%> 


2*  d?  Rtf 

■  I  T  l 

♦f *0  0j“6i  R*Rj 


1-j 


where  P  ■  1  R  >  cutoff  pertaining  to  direction  6,* 
n0R<"  "  •  •  • 


with  Rc(0i*^j)  the  cutoff  pertaining  to  the  direction  0n 

refers  to  the  zenith  angle  range  lying  within  the  ri^1  sene  (extending  from 
$1  to  62) .  The  function  C  is  the  sane  as  that  in  equation  2.1,  now 
integrated  over  all  azimuths.  If  desired,  limited  ranges  of  the  parameter 
4  can  be  introduced,  in  which  case  an  azimuth  dependence  exists. 


The  cutoff  distribution  functions  calculated  far  the  sites  listed  in 
Table  2.2  are  presented  in  Figures  2.1-8.  The  zenith  angle  and  latitude 
dependence  of  the  functions  are  well  displayed  in  these  figures,  which  are 
arranged  in  order  of  decreasing  magnetic  latitude. 

It  is  believed  that  these  functions  have  basically  the  correct  form, 
in  spite  of  the  simplified  approach  to  the  calculations.  It  is  worth 
noting  that  any  step  to  improve  the  precision  of  the  functions,  either  by 
introducing  real  field  cutoffs,  or  by  talcing  into  account  the  transverse 
momentum  in  the  primary  ooamic  ray  interactions,  would  require  an  amount  of 
oenpuber  time  greater  by  many  orders  of  magnitude.  That  that  neutrino  event 
rates  observed  by  the  particle  detectors  are  extremaly  low,  in  any  case, 
argues  far  the  simpler  approach  employed  at  this  time. 


The  effect  of  the  simplifying  assumptions  an  the  functions  may  be 
anticipated.  Note,  far  example,  that  for  zones  1  through  5,  the  median 
value  of  the  functions  lie  at  a  rigidity  value  corresponding  to  the 
vertical  Sterner  cutoff  value.  This  fact  suggests  that,  as  a  first  order 
correction  to  take  into  aooount  the  disparity  between  the  Starmsr  end  real 
field  cutoffs,  these  curves  could  be  displaced  to  positions  for  which  the 
median  value  of  the  functions  lie  at  the  real  field  vertical  effective 
cutoff  pertaining  to  each  individual  site. 
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If  a  more  realistic,  three  dimensional,  neutrino  production  model  were 
to  toe  utilized,  a  "smearing"  of  the  functions  would  be  expected,  which 
would  probably  result  in  less  well  defined  zenith  angle  dependence,  on  the 
other  hand,  the  effect  of  the  Earth's  shadow,  and  of  secondary  particle 
deflection  over  large  atmospheric  path  lengths,  could  be  expected  to  cause 
the  functions  applying  to  near  horizontal  directions  to  extend  to  rigidity 
values  appreciably  greater  than  the  curves  indicate. 

Gaisser  et  al. ,  1983,  have  recently  used  the  cutoff  distribution 
functions  (as  ocrputed  in  a  form  relating  to  upwards  and  downward  directed 
detector  response  cones  of  various  half  angles)  in  calculating  the  flux  of 
atmospheric  neutrinos.  Ihey  have  shewn  that  the  geomagnetic  effect  has  a 
very  significant  effect  on  the  expected  up-down  ratios  of  neutrinos  of  the 
two  flavors,  and  hence  on  the  interpretation  of  the  experimental  data  from 
the  large  detectors. 

TABLE  2.2 

locations  for  which  cutoff  distribution  functions  have  been 
calculated.  A  single  calculation  has  been  carried  out  for  the  two  tunnel 
sites  on  the  French/I talian  border  because  the  two  experiment  locations  are 
sufficiently  close  spaced  to  possess  essentially  the  same  cutoff 
distribution  functions. 


RlG'DITY  (GV) 

FIGURES  2. 2-2. 5  (in  descending  cedar) 

Integral  cutoff  distribution  functions  for  the  locations  shown  in 
table  2.2.  For  aach  location  a  sat  of  curves  is  present  ad,  each  wwhsr  of 
which  corresponds  to  a  different  solid  angle  sons  and  represents  tbs 
percentage  of  the  total  solid  angle  of  the  cone  accessible  to  neutrinos 
that  are  descended  from  charged  primaries  of  any  given  rigidity.  3he 
zenith-angle  extant  of  each  zona  is  identified  in  table  2.1. 
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0.1 


1  10 
RIGIDITY  (GV) 


FIGURES  2. 6-2. 9  (in  descending  order) 
See  caption  for  preceding  figures. 
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CHAFER  3:  THE  EFFECT  OF  MAGNETIC  CUTOFFS  CM  THE  FIDX  OF  CHARGES 

PRIMARY  COSMIC  HAYS  WITHIN  THE  MAIN  FIELD  OF  JUPITER 


3.1  INIHOCUCnCN 

With  the  prospect  of  the  Galileo  spacecraft  eocn  to  be  placed  in  orbit 
around  Jupiter,  there  is  a  direct  interest  in  estimating  the  charged 
particle  fluxes  in  the  orbital  snvircnsant.  Although  Stonier  cutoffs  could 
have  been  used  to  approximate  cutoff  values  (as  dons  for  Jupiter's 
magnetosphere  by  Cooke,  1974)  a  higher  precisian  was  desired  at  the  present 
time.  In  the  present  investigation,  therefore,  "real"  main  cutoff  values 
ware  derived  using  visual  inspection  of  trajectories  plotter  by  ocnpubsr. 

In  particular,  trajectories  were  presented  cn  a  video  screen  far  a  range  of 
rigidity  values,  corresponding  to  particular  locations  and  directions  of 
arrival  in  JUpiter' s  field.  The  farm  of  the  trajectory  corresponding  to 
each  set  of  conditions  was  inspected  as  the  rigidity  value  progressively 
refined  until  the  simplest  possible  quasi-bound  periodic  orbit  was  found  to 
be  associated  with  the  approach  trajectory.  The  rigidity  value  thus  found 
corresponds  to  the  main  cutoff  value  (this  association  of  the  simplest 
quasi-periodic  orbit  with  the  main  cone  is  discussed  further  in  chapter  6, 
in  relation  to  the  geomagnetic  field) . 

Techniques  described  in  the  preceding  chapter  have  been  applied  to 
calculate  the  effect  of  the  magnetic  cutoffs  cn  the  primaries  within  the 
magnetic  field  of  Jupiter,  expressed  by  means  of  cutoff  distribution 
functions  (CEF's) .  Cutoff  distribution  functions  have  been  derived  for 
primary  cosmic  rays  having  aoosss  to  a  4.5  Jupiter  radius  satellite  orbit, 
far  longitudes  over  the  latitude  range  +20  to  -2CP.  These  functions 
represent  the  fraction  of  tha  total  4*  steradian  solid  angle  which  is 
accessible  to  charged  primaries  of  any  given  rigidity. 


3.2  DISCUSSION 

Let  tha  differential  intensity  of  charged  primaries  reaching  a  qivm 
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point,  in  the  presence  of  Jupiter's  field  be  A(R,x,4,r)  (where  R  is 
rigidity,  x  and  <j>  are  the  latitude  and  longitude,  and  r  is  the  radial 
position  of  the  point  in  the  field).  If  B(R,x,4,r)  is  the  intensity  of 
charged  primaries  that  would  exist  in  the  absence  of  the  field,  then  the 
functions  A  and  B  can  be  related  as  follows: 

A(R,X,^,r)  ■  B(R,A,*,r)  C(R,x,4,r) 

where  C(R,x,<f>,r) ,  the 

cutoff  distribution  function,  expresses  the  effect  of  the  magnetic  cutoffs. 
It  defines,  for  a  given  location,  the  fraction  of  the  total  solid  angle 
accessible  to  charged  primaries  of  any  given  rigidity. 

The  cccputer-based  procedure  described  in  chapter  2,  there  enplcyed 
far  calculations  in  the  geomagnetic  field,  is  used  to  calculate  these 
functions,  using  a  pre-calculated  set  of  main  cutoff  values.  Ihese 
cutoffs,  which  were  calculated  for  a  spaced  set  of  zenith  and  azimuth 
angles  pertaining  to  a  grid  of  latitude  and  longitude  values,  specify,  far 
any  given  direction  and  location,  the  lower  limit  of  full  accessibility  of 
primaries  from  outside  the  field.  As  briefly  discussed  earlier,  this  lower 
limit,  the  "main  cutoff"  (see  chapter  1) ,  was  identified,  in  each  case,  by 
locating  the  rigidity  value  far  which  the  simplest  bound  periodic  orbit  was 
able  to  exist.  Uve  O4  mathematical  representation  of  JUpiter' e  main  field 
(Acuna  and  Ness,  1976)  was  used  in  making  the  cutoff  calculations.  During 
these  oonputaticns  it  was  observed  that  relatively  few  allowed  trajectories 
were  encountered  at  rigidities  below  the  main  cutoff,  tending  to  validate 
the  use  of  the  main  cutoff  to  represent  the  lower  limit  of  accessibility. 

At  4.5  Jtpiter  radii,  to  which  the  cutoff  values  pertain,  the  effect  of  the 


FIGURE  3.1 

General  fora  of  the  cutoff 
distribution  function.  The  rigidity 
values  R^,  and  fy,  used  to 
characterize  functions  of  this  kind, 
are  shown. 
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"shadow**  of  the  solid  planet  is  restricted  to  a  very  snail  (-1%)  effect  at 
very  high  rigidities. 

In  detail,  the  derivation  of  the  CDF's  for  any  location  involves 
summation  of  the  effects  of  the  cutoffs  over  many  small  solid  angle 
elements,  using  cutoff  values  interpolated  from  title  stored  set.  Because 
charged  particle  aooess  varies  from  location  to  location,  the  resulting  CDF 
is  different  far  each  different  location.  The  representation  of  the  large 
number  of  Individual  functions  pertaining  to  a  wide  range  of  latitude  and 
longitude  values  is  clearly  impractical  here,  so  an  alternate  graphical 
representation  has  been  employed,  which  uses  three  contour  plots  to  show 
the  variation  in  the  value  of  each  of  three  parameters  which  ere  here  used 
to  characterise  the  function.  Figure  3.1  shows  the  typical  form  of  e  CDF. 
Rl,  %,  and  %  are  three  rigidity  values  which  respectively  represent  the 
lower  limit,  the  rigidity  at  which  the  function  has  value  50%,  and  the 
upper  rigidity  limit  of  the  function. 


3.3  RESULTS: 


Figures  3.2A,  3.2B,  and  3.2C  show  contour  plots  of  the  values  of  Rl, 
%,  and  F(j,  respectively,  far  all  longitudes  over  the  latitude  range  +20  to 
>20°.  At  any  desired  location  the  values  of  the  three  parameters  can  be 
deduced,  and  used  to  construct  the  appropriate  CDF.  The  contour  plots 
pertain  to  a  4,5  Jupiter  radius  shell.  Approximate  functions  far  other 
radius  values  may  be  deduced  by  invoking  an  inverse  square  law  of  variation 
of  rigidity  value. 

Finally,  the  flux  at  any  given  location  can  be  deduced  by  folding  the 
appropriate  CDF  with  the  primary  spectna.  If  it  is  desired  that  energy  be 
used  rather  than  rigidity  as  the  spectral  variable,  then  at  the  rigidity 
values  involved  the  conversion  Is  siaple.  Proton  energy  is  almost  exactly 
numerically  equal  to  rigidity;  alpha  nucleon  energy  is  numerically  equal  to 
half  the  rigidity  value. 
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FIGURE  3.2A  (top  panel) 

Oontour  nap  of  the  value  of  the  parameter  %,  one  of  the  three 
rigidity  values  used  to  characterize  the  cutoff  distribution  functions. 

FIGURE  3.2B  (central  panel) 

Contour  nap  of  the  value  of  Jfy. 


FIGURE  3.2C  (lowest  panel) 


Oontour  nap  of  the  value  of 
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CHAPTER  4:  TRAJECTORY  PARAMETERIZATION 


4.1  INIKODUCTXCN 

Prior  to  the  development  of  the  techniques  presented  in  this  report 
the  stixJy  of  oosmic  ray  access  into  planetary  magnetospheres  (mainly  that 
of  the  Earth)  had  been  carried  cut  by  methods  which  could  broadly  be 
divided  into  two  types.  The  earlier,  analytic  or  quasi-analytic  approaches 
(including  the  work  of  Stonier,  1930;  Lemaitre  and  Vallarta,  1936a, b;  for 
exanple)  were  limited  to  consideration  of  field  models  which  possessed 
axial  symnetry.  These  approaches,  whilst  significantly  limited  in  the 
extent  to  which  they  could  represent  the  real  geomagnetic  field, 
nevertheless  yielded  very  valuable  insights  into  the  phenomenology  of 
oosmic  ray  access  to  points  within  the  field. 

The  later  work,  utilizing  the  speed  and  power  of  the  digital  computer, 
involved  direct  testing  of  access  by  systematic  treeing  of  trajectories  at 
close  spaced  rigidities  far  the  :  ange  of  final  directions  at  the  point  of 
interest  within  the  magnetosphere.  This  technique  is  capable  of  producing 
cutoffs  to  as  high  a  precisian  as  required  (limited  cnly  by  the  accuracy  of 
the  modelling  of  the  field,  atmosphere  etc. ) ,  but  is  not  capable  of  giving 
ary  insight  into  the  phenomenology  of  oosmic  ray  access,  and  has  proved  to 
be  a  very  inefficient  means  of  mapping  cutoffs.  In  practice,  trajectories 
have  to  be  traced  for  every  possible  angle,  position  and  rigidity  in  any 
location,  in  order  to  build  up  a  picture  of  cosmic  ray  access  in  the  given 
situation. 

The  research  summarized  in  this  report  has  involved  developing 
techniques  far  numerically  characterizing  trajectories  in  order  to  allcw  a 
oonputer  to  relate  trajectories  over  any  ccntinuun  of  rigidity,  angle,  and 
position,  far  exzmple.  This  technique,  referred  to  as  "Trajectory 
parameterization",  allows  full  insist  into  the  relationship  betwesn 
trajectory  characteristics,  such  as  local  points  of  maximum  and  minima 
altitude,  equator  crossings,  and  so  on,  and  the  types  of  aaosss  regions 
sumnarized  in  the  Chapter  1. 
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The  development  of  the  techniques  passed  through  a  number  of  different 
phases,  before  reaching  the  fully  developed  form  finally  achieved.  This 
chapter  reviews  this  process  of  development. 

4.2  RELATICN  BETWEEN  TRAJECTORY  CHARACTERISTICS  AND  PENUMERAL  FEATURES 

The  iterative  digital  ocnputer-based  technique  determines,  by  direct 
testing,  the  ability  of  particles  of  any  particular  rigidity  to  arrive  at  a 
given  site  in  &  given  direction.  A  "negative"  trajectory  trace  is  carried 
out,  with  the  computation  starting  at  the  site,  in  the  direction  of 
interest.  Failure  of  the  trajectory,  as  calculated  in  the  reverse 
direction,  to  travel  successfully  to  the  outer  boundary  of  the  geomagnetic 
field  indicates  that  a  cosmic  ray  particle  would  be  unable  to  arrive  at  the 
site  under  the  given  conditions. 

For  any  particular  direction  at  a  given  site  there  is  generally  a 
complex  intermingling  of  escaping  (in  the  negative  time  sense) ,  non¬ 
escaping  (including  those  which  intersect  the  surface  of  the  planet) 
trajectories,  which  gives  rise  to  a  complicated  pattern  of  allowed  and 
forbidden  rigidity  values  (i.e.  the  penumbra) .  The  complexity  arises  out 
of  the  characteristic  trajectory  deflection  produced  by  the  Lorentz  farce 
which  affects  the  charged  particle  as  it  moves  through  the  magnetic  field. 
Calculations  shew  that  particles  destined  to  arrive  at  any  given  site  on 
the  Earth  with  rigidity  values  within  the  penumbra  typically  travel  from 
east  to  west  towards  the  site,  and  swing  back  and  forth  across  the  magnetic 
equator.  In  the  longitude  range  in  the  proximity  of  the  arrival  site  the 
trajectories  become  tied  to  the  local  field  line  bundle,  tending  to  loop 
within  this  bundle  in  the  final  stages  of  their  flight  towards  the  site. 

Figure  4.1  shows  a  trajectory  in  which  this  behavior  is  very  clearly 
displayed.  The  trajectory  (as  traoed  "negatively**  -  i.e  away  from  the 
point  of  arrival)  moves  from  west  to  east,  and  swings  repeatedly  from  one 
side  of  the  geomagnetic  equator  to  the  other,  with  the  altitude  changing 
continuously.  If,  in  one  of  the  low  points,  the  surface  of  the  Earth  is 
encountered,  cosmic  ray  access  via  that  trajectory  would  not  be  allowed. 


w 

OS 


H 

X 

o 

HH 

W 

X 


at 

» 


3  S 


'  >2  *2. 


H  IN.  »V#  IK.  lift-  >«».  m  MO. 


Latitude 

40.00 

Longitude 

180.00 

Zenith 

100.0 

Azimuth 

140.0 

Altitude 

400 

Re-En t  Alt 

30 

Rigidity 

4.849 

Path  Length 

34.666 

NStep 

2328 

Min  Alt  (KM) 

23.6 

NMax 

29 

Re-Entrant 

Latitude 

43.48 

Longitude 

504.40 

ir  «r 

I; 


IN  l%0  IK  fit*  !«•  m 


K  !•' 


»* 


»r 


LONGITUDE  RELATIVE  TO  SITE  (DEG) 


HEIGHT  (RE) 


FIGURE  4.1 

A  trajectory  displaying  narked  quasi-bound  periodic  orbit  behavior. 
Note  that  the  plots  use  geographic  latitude  and  longitude  scales.  The 
position  of  the  cosmic  ray  equator  nay  clearly  be  seen  from  examination  of 
the  distortion  of  the  gross  trajectory  fora  in  the  latitude  vs.  longitude 
plot. 
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The  intersections  demonstrably  responsible  for  the  failure  of 
particles  to  enter  along  a  particular  trajectory  may  in  principle  be 
classified  according  to  the  section  of  the  trajectory  in  which  they  occur. 
Consider  such  trajectories  as  being  followed  in  the  "negative"  direction  - 
traced  outward  away  from  the  site.  The  trajectories  are  seen  to  pass 
through  loops  and  other  contortions  which  take  the  "particle"  alternately 
towards  and  away  from  the  planet.  If,  in  one  of  the  approaching  sections, 
the  planet's  surface  (i.e.  the  top  of  the  atmosphere,  in  the  case  of  the 
Earth)  is  intersected,  then  that  particular  trajectory  will  be  forbidden  to 
particles  of  the  given  rigidity  from  outside  the  magnetic  field. 

The  point  along  the  trajectory  at  which  reentry  occurs  characterizes 
the  allowed/ forbidden  structure  of  the  penumbra  in  the  physical 
locaticn/direction/rigidity  space  of  the  particle.  Very  clearly  defined 
regions  of  this  space  may  be  distinguished.  First  of  all,  we  will  consider 
those  produced  by  intersection  with  the  planet  of  the  particle  trajectories 
within  the  local  field  line  bundle.  The  traditional  shadow  cone  forbidden 
region  is  due  to  intersection  of  trajectories  in  the  immediate  locality  of 
the  site  in  question  before  even  a  single  loop  is  performed.  Other  shadow 
cone  regions  are  associated  with  the  intersection  of  the  second,  third,  or 
etc.  loop  with  the  surface  of  the  planet. 

Intersection  of  these  loops  with  the  surface  of  the  planet  occurs  in 
the  general  vicinity  of  the  site.  Because  of  the  short  path  traversed  by 
the  ocnputatianal  "particle"  in  traveling  from  the  site  to  the  point  of 
intersection  of  the  loops  with  the  surface  of  the  planet,  and  because  the 
paths  are  restricted  to  relatively  low  altitudes,  the  pemmibral  structures 
so  produoed  are  relatively  stable  and  field  independent. 

On  the  other  hand,  if  the  negative  "particle"  travels  successfully 
through  the  local  field,  negotiating  the  loops  (if  present)  without 
intersecting  the  planet's  surface,  then  its  motion  takes  on  a  more  broadly 
identifiable  character.  Such  trajectories  are  seen  to  travel  back  and 
forth  between  the  northern  and  southern  magnetic  hemispheres,  while  moving 
from  west  to  east.  This  portion  of  the  trajectories  will  be  referred  to  as 
the  periodic  orbit  part.  The  trajectory  may  possess  zero,  one,  or  wars 
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loops  in  each  section  (we  define  a  section  as  ary  individual  part  of  a 
trajectory  lying  to  one  side  or  other  of  the  magnetic  equator) ,  depending 
an  the  nature  of  the  particular  trajectory*  Within  the  loops  the  charged 
particle  is  able  to  approach  the  planet  most  closely,  and  the  trajectory 
may  intersect  the  planet's  surface  at  one  of  these  points.  Although 
trajectories  in  the  immediate  physical  locaticrVdirecticrv'rigidity  apace 
will  also  be  affected,  the  fact  that  the  path  between  the  site  and  the 
reentry  point  is  relatively  lcng  means  that  the  forbidden  space  tends  to  be 
of  restricted  extent,  and  to  be  parameter  and  field  sensitive. 

Because  of  the  greatly  increased  number  of  near  approach  points  in 
long  multiple  section  trajectories,  and  the  greater  trajectory  path  lengths 
involved,  the  panunfaral  structure  so  produced  is  very  much  more  finely 
structured,  oanplax  and  field  sensitive  than  the  shadow  cone  structures. 
Nevertheless,  it  is  in  principle  possible  to  classify  the  penunbral 
structure  detail  according  to  the  section  of  the  trajectory  in  which  the 
reentry  has  occurred,  and  to  explore  the  sensitivity  of  the  structure  to 
the  field  and  parameter  changes.  Systematic  labelling  of  the  trajectory 
structures  would  open  the  way  for  investigating  the  nature  of  the  penuntora, 
and  the  other  allowecVforbidden  structures  (as  reviewed  in  Chapter  1) . 

As  a  first  step  towards  implementing  such  a  scheme  far  use  in  an 
automated  computer-based  system,  trajectories  were  characterized  in  terms 
of  the  number  of  loops  performed  in  each  of  the  trajectory  sections.  A 
standard  trajectory  tracing  program  was  modified  to  allow  the  appropriate 
detection  and  counting  of  trajectory  loops. 

that  is  a  loop?  A  loop  is  considered  to  have  been  executed  whan  the 
curvature  vector  associated  with  the  trajectory  has  rotated  through  360 P. 
Figure  4.2  shows  a  reentrant  trajectory  within  the  geomagnetic  field  which 
contains  six  loops,  three  very  obvious  ones  (seen  clearly  in  the  altitude 
vs.  longitude  plot) ,  and  three  othare  -  cane  near  the  site,  one  more  just 
before  the  first  crossing  of  the  geomagnetic  equator,  and  another  in  the 
last  section  of  the  trajectory,  in  its  passage  toeerds  the  mirfaoe  of  the 
earth.  In  practice  loops  such  as  these  are  detected  automatically  ly 
ocpputar  during  the  trajectory  trace.  The  number  of  loops  in  each  section 
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LATITUDE  (DEG)  «"  HEIGHT  (RE) 

6C  *  -30.  0.  30.  60.  SO-  1.0  '  10'  10* 


FIGURE  4.2 

A  sample  trajectory,  showing  six  loops  (categorized  by  ocnputer  as  a 
2211  reentrant  trajectory) . 
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is  counted  as  the  trace  proceeds,  and  is  stored  for  printing  at  the  and  of 
the  calculation.  The  order  of  the  loop  numbers  as  they  axe  listed  (in 
figure  4.3,  far  example)  is  such  that  the  first  number  refers  to  the 
section  closest  to  the  site. 

Figure  4.3  shows  a  plot  based  on  the  trajectory  loop  numbers.  This 
plot  represents  a  section  of  the  terrestrial  penumbra  far  an  azimuth  range 
of  129-145°  geographic  for  the  given  situation  (100°  zenith  angle,  400  lan 
altitude,  40°  latitude,  180°  longitude,  field  model  IGRF  1975.0  modified  to 
Epoch  1980.0) .  Trajectory  loop  numbers  are  presented  for  1%  rigidity 
intervals  in  the  range  10.025  to  4.054  GV.  This  situation  was  investigated 
as  part  of  another  study,  and  is  not  represented  as  having  any  particular 
significance  to  this  present  discussion.  Although  the  particular  structure 
(including  the  "overhung*'  shadow  ocne,  made  up  of  unlooped  reentrant 
trajectories)  is  in  detail  characteristic  only  of  the  given  site  and 
conditions,  the  applicability  of  the  category  scheme  to  a  cccplex  situation 
such  as  this  is  well  illustrated. 

Inspection  of  the  listed  categories  shows  that  there  exists  a  well- 
ordered  structure  of  penunibral  bands,  even  in  the  continuum  of  forbidden 
trajectories  below  the  few  allowed  trajectories  in  the  region  of  8-9  GV. 
Bands  have  been  shaded  to  emphasize  the  band  structure.  Calculated  loop 
numbers  fox  intermediate  azimuths  (not  shown)  have  been  vised  to  help 
resolve  the  complexity  in  sane  areas.  At  the  It  rigidity  interval  utilized 
here  the  wider  bands  attributable  to  short  and  medium  range  reentrant 
trajectories  (involving  trajectories  having  up  to  approximately  3 
equatorial  crossings)  show  clearly.  With  decreasing  rigidity  an  orderly 
succession  of  bands  appears,  produced  by  trajectories  possessing  greater 
and  greater  number  of  loops.  At  any  point  the  low  order  bands  due  to  short 
range  earth  intersections  overlie  the  higher  order  bands  resulting  from 
intersections  further  along  the  trajectories.  An  apparent  continuum  of 
forbidden  trajectories  may  in  fact  consist  of  many  overlapping  bands  (where 
the  latter  may  be  infinitesimally  narrow) .  calculations  carried  cut  for 
progressively  finer  rigidity  intervals  (and  for  finer  azimuth  spacing)  tend 
to  reveal  progressively  greater  detail  in  the  high  order  penumbral  band 
structure. 
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FIGURE  4.3 

Listed  ccrputer  derived  loop  numbers  (truncated  to  a  maximum  of  12 
digits)  for  a  sample  penuntoral  region  at  100°  zenith  angle  (the  location 
etc. ,  is  defined  the  text) .  Clearly  distinguished  penuntoral  bands  are 
numerically  identified  and  shaded  (heaviest  shading:  1  trajectory  equator 
crossing,  lightest:  4  crossings) .  The  asterisks  identify  reentrant  orbits 
Unshaded  areas  in  the  penumbra  may  contain  either  forbidden  bards  or 
allowed  trajectories. 


This  technique  revealed  seme  very  interesting  insights  into  sene 
particular  aspects  of  the  main  oane  and  penumbra. 

a)  In  the  calculated  penumbra  it  is  the  forbidden  bands,  not  the 
allowed  bands,  that  are  defined.  An  allowed  band  can  only  exist 
when  no  obscuring  forbidden  bands  lie  at  the  rigidity  in  question 
(analogous  to  looking  at  a  distant  view  between  the  trunks  of  a 
forest  of  trees  -  one  sees  the  view  when  there  are  no  trunks  in 
the  way ;  it  only  takes  one  trunk  to  obscure  the  view,  no  natter 
hew  many  gaps  there  are  also  in  the  line-of-sight) . 

b)  Penumbral  structure  produced  by  earth  intersections  in  sections 
close  to  the  site  are  much  wider  and  more  stable  than  those 
resulting  firm  longer  range  reentrant  trajectories. 

c)  The  offset  of  the  basic  geomagnetic  dipole  manifests  itself  In 
differences  in  the  degree  of  overlapping  and  nixing  of  penumbral 
bands  at  different  points  on  the  Earth's  surface. 

d)  Higher  order  penumbral  bawds  (those  possessing  several  sections) 
have  greater  sensitivity  to  geomagnetic  field  perturbations  and 
to  changes  in  field  node!  parameters  than  do  those  of  lmer 
order,  and  so  are  relatively  unstable. 

There  is  a  very  high  degree  of  order  in  the  penumbra  which  is  usefully 
displayed  in  diagrams  of  the  type  that  figure  4.3  represents.  This  kind  of 
plot  has  to  be  prepared  by  inspection  of  the  arrays  of  loop  numbers,  a  very 
inefficient  process.  Another  kind  of  plot  can  better  be  used  to  chew  the 
penumbral  order  -  one  involving  a  plot  of  trajectory  feature  positions, 
which  is  prepared  completely  by  oenputar. 


4.3  SUMMARY  PLOT  REWESENIATICN  OF  THE  PENUMBRA 

The  "Summary  Plot"  representation  is  nose  efficient,  and  a  somewhat 
more  refined  way  of  representing  trajectory  aooess  information. 
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Trajectories  are  effectively  summarized  in  terms  of  plots  of  the  positions 
of  the  trajectory  "cardinal  points'*.  These  cardinal  points  are  loops#  lew 
points  and  magnetic  equator  crossing  points.  These  points  change  position 
(both  path  length  and  altitude)  in  a  systematic  manner  with  changing 
rigidity.  The  cardinal  point  positions  may  be  plotted  to  produce  the  so- 
called  summary  plot,  in  which  for  each  of  the  types  of  cardinal  points  the 
point  positions  are  presented  for  a  range  of  values  of  rigidity  (or  for  a 
range  of  whatever  parameter  is  of  particular  interest) .  These  plots  are 
prepared  on  a  printer,  using  printed  characters  positioned  on  a  page  so 
that  their  position  is  related  to  the  path  length  along  the  trajectory,  as 
traced  in  the  "negative  time"  direction,  away  from  the  final  arrival 
point) .  The  character  used  to  represent  each  point  represents  useful 
information  about  the  point;  for  example,  altitude  in  the  case  of  low 
points,  and  "loop  development"  in  the  case  of  the  loop  position  plot. 

Figure  4.4  shews  a  typical  summary  plot  set  (this  is  a  direct 
reproduction  of  a  microfiche,  as  reproduced  from  Ooctke  et  al. ,  1981). 

Cooke  and  Bredesen,  1981,  presented  another  set  of  summary  plots  which, 
like  these  reveal  a  very  highly  ordered  structure.  In  figure  4.4  the 
penumbral  forbidden  bands  are  clearly  shown  as  points  where  the  low  points 
intersect  the  reference  height  corresponding  to  what  is  deemed  to  be  the 
grazing  height  in  the  atmosphere.  The  penumbral  bands  are  directly 
displayed  on  the  plot  which  shows  the  position  of  the  intersection  points. 
Note  that  the  edge  of  the  Stormer  cane  is  displayed  clearly  in  this  plot, 
as  the  boundary  belcw  which  all  trajectories  display  bound  periodic 
behavior.  After  the  onset  of  this  behavior,  the  trajectories  are 
carpi etely  bound,  and  cannot  escape  the  field  (and  therefore  cosmic  rays 
outside  the  field  with  the  particular  rigidities  are  unable  to  gain  access 
to  the  particular  point  and  direction  involved) . 

The  technique  is  employed  of  continuing  the  trajectory  traces  beyond 
these  intersection  points  so  that  structures  may  be  traced  from  one  side  to 
the  other  of  the  range  of  rigidity  in  which  a  forbidden  band  lies.  It 
should  be  clearly  understood,  however,  that  the  intersection  trajectories 
do  not  represent  "read"  trajectories,  tut  are  merely  a  useful  aid  to  the 
tracing  of  continuous  structures. 
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FIGURE  4.4 

Suimnary  plot  set  of  Cooke  and  Bredesen,  1981,  representing  positions 
of  equatorial  crossings,  loops,  low  points,  and  eartVatmosphere 
intersections  (paired  columns  from  left  to  right) ,  for  a  range  of  rigidity 
values  from  15.00  GV  down  to  5.00  GV  in  0.05  GV  increments.  These 
calculations  were  carried  cut  for  vertical  incidence  at  the  location  with 
geographic  latitude  -10°  and  longitude  270°,  and  400  km  altitude?  far 
magnetic  field  IGRF  1975.0  updated  to  1980.0.  The  paired  diagrams  show 
detail  for  trajectory  path  lengths  of  0-12  and  12-24  earth  radii  The 
distinct  change  in  behavior  of  the  longer  trajectory  sections,  which  occurs 
at  7.20  GV,  is  interpreted  as  being  produced  by  transit  of  the  Stonier  cere 
edge. 
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The  scale  of  the  trajectory  cardinal  point  structures  and  the 
penumbral  structures  varies  ever  a  vide  range  -  those  a  short  range 
trajectory  sections  say  exist  for  many  GV,  while  at  long  trajectory  ranges 
the  structures  say  exist  only  fleetingly  and  say  be  represented  by  only  a 
single  point  an  a  summary  plot,  if  seen  at  all.  Nevertheless,  when 
examined  at  appropriately  fine  rigidity  intervals  they  show  the  same 
characteristic  forms  as  the  larger  scale,  short  range,  structures. 


4.4  AUTOMATED  TRAJECTORY  PARAMETERIZATION 

As  a  final  stage  in  the  use  of  trajectory  parameters  to  aid  in  the 
exploration  of  cosmic  ray  access  to  points  within  the  geomagnetic  field, 
strategies  have  been  developed  and  incorporated  into  computer  programs  that 
locate  the  trajectory  cardinal  points  and  then  correlate  the  position  of 
these  points  with  change  in  ary  of  the  parameters  involved  in  the 
computations,  in  order  to  map  the  features  associated  with  the  rard-lnui 
points  with  any  required  range  of  parameter  value.  The  new  approach,  while 
utilizing  the  speed,  efficiency  arid  "real"  geomagnetic  field  modeling 
capabilities  of  the  digital  computer,  yields  an  analytical  insight 
equivalent  to  the  earlier  and  elegant  approaches  of  Stormer  and  Lemaitre 
and  Vallarta.  The  technique  is  applicable  to  cosmic  ray  aooess  into  any 
planetary  magnetosphere,  although  the  discussion  here  is  presented  in  terms 
of  the  situation  in  the  geomagnetic  field,  in  relation  to  which  the 
technique  was  developed. 

Like  the  basic  digital  computer  iterative  penumbral  mapping  procedure 
out  of  which  it  has  grc*m  (Shea  et  al. ,  1965),  the  parameterization 
approach  involves  the  tracing  of  trajectories  in  the  reverse,  "negative 
time",  direction,  out  away  from  the  point  in  the  magnetic  field  for  which 
the  penumbra  is  to  be  studied.  Similarly,  the  input  parameters  required 
for  the  trajectory  traces  include  those  identifying  the  location  and 
direction  of  interest  at  the  site  -  latitude,  longitude,  altitude,  zenith 
and  azimuth  -  and  other  parameters  such  as  particle  rigidity,  grazing 
altitude  (altitude  above  the  Earth's  surfaoe  at  which  a  trajectory  is 
deemed  to  graze  the  atmosphere) ,  magnetic  field  parameters,  and  time. 


61 


Whereas  the  standard  penumbral  mapping  technique  involves  a  systematic 
testing  of  whether  individual  trajectories  are  allowed  or  forbidden  for 
sets  of  discrete  parameter  values  (typically  for  varying  rigidity) ,  the 
"parameterization"  approach  takes  advantage  of  the  fact  that  trajectories 
possess  a  continuity  of  form  over  variation  in  any  parameter  in  the  set  to 
efficiently  trace  particular  penumbral  features,  or  to  elucidate  their 
characteristics. 


Although  the  trajectory  configuration  may  in  same  situations  alter 
exceedingly  rapidly  with  variation  of  a  given  parameter,  the  changes  in 
form  take  place  in  a  continuous  manner.  Hie  trajectory  "cardinal  points" 
are  found  to  change  position  in  a  quite  systematic  manner  with  changes  In 
any  parameter.  Figure  4.5,  an  enlarged  section  of  one  of  the  trajectory 
summary  plots  of  Cooke  and  Bredesen,  1981,  illustrates  this  clearly. 

As  recognized  by  Lemaitre  and  Vallarta  (1936a) ,  the  generators  of 
penumbral  structures  (i.e. ,  the  trajectories  defining  the  edge  of  the 
allowed  regions,  lying  at  the  transition  from  allowed  to  forbidden  regions) 
are  allowed  trajectories  that  graze  the  planet  (Earth's  atmosphere)  at  sane 
point  along  their  length. 

Assume  that  the  parameters  belonging  to  a  transition  trajectory  in  a 
given  situation  have  been  established  -  the  trajectory  has  been  traced  and 
found  to  graze  the  planet  at  a  low  point  associated  with  one  of  the  loops 
along  the  trajectory.  Varying  the  value  of  any  one  parameter  over  a  small 
interval  will  tend  to  alter  the  lew  point  altitude  so  that  the  trajectory, 
initially  grazing,  is  new  clearly  allowed  or  forbidden  (depending  on 
whether  the  low  point  rises  or  falls) .  Normally  it  is  possible  to  vary  the 
value  of  one  of  the  other  parameters  in  order  to  again  make  the  same 
trajectory  low  point  graze  the  surface  even  though,  of  ccurse,  the  overall 
trajectory  form  will  have  changed  slightly  with  the  matched  change  in  the 
two  parameters. 


Biis  operation  of  matching  the  change  in  two  parameters  in  order  to 
maintain  trajectory  grazing  is  the  key  to  the  parameterization  technique. 

A  family  of  techniques  can  be  implemented  through  the  suitable  use  of  this 
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basic  step  -  mapping,  determination  of  the  stability  of  penumbral 
structures,  •♦migration"  of  structures  over  ranges  of  parameters  (in  order 
to  establish,  for  example,  change  in  a  structure  along  the  path  of  a 
satellite  in  orbit) ,  to  name  a  few  of  the  techniques  that  have  already  been 
implemented  in  the  computer  program  developed  for  use  in  the  study  of  the 
penumbra. 


The  parameterization  approach  depends  critically  an  the  ability  of  the 
computer  program  to  reliably  identify  particular  cardinal  points  and  to 
follow  these  in  the  presence  of  continued  change  in  trajectory 
configuration,  even  though  the  trajectories  are  generally  complex  and 
possess  a  changing  number  of  low  points  and  loops.  The  key  to  the 
automatic  tracking  of  trajectory  cardinal  points  is  the  vise  of  position  and 
"development"  parameters.  In  practice  the  position  of  the  given  cardinal 
feature  (normally  a  low  point)  along  a  trajectory,  as  measured  from  the 
trajectory  starting  point,  is  monitored,  and  the  observed  changes  in 
position  are  related  to  the  change  in  the  variable  input  parameter. 
Interpolation  or  extrapolation  is  then  employed  to  calculate  where  the 
feature  would  be  expected  to  lie  under  the  modified  set  of  conditions.  That 
a  certain  feature  lies  close  to  the  anticipated  position  is,  however,  not  a 
sufficient  condition  for  unambiguous  identification  of  the  trajectory 
features  of  particular  interest.  The  "development  parameter",  related  to 
the  degree  of  development  of  each  feature  (and  in  part  related  to  the 
trajectory  curvature  at  the  point),  is  computed  for  each  feature  along  a 
trajectory.  This  parameter,  while  varying  relatively  slowly  with  change  in 
the  input  parameter,  differs  significantly  from  feature  to  feature,  and  so 
can  be  used  as  an  identification  label  to  aid  recognition  of  particular 
features.  The  application  and  use  of  both  the  position  and  development 
parameters  has  been  found  sufficient  to  allow  the  computer  program  to 
"recognize"  features,  and  to  monitor  their  role  in  contributing  to 
penumbral  structures. 

Because  the  iterative  trajectory  trace  procedure  is  Inherently  a 
noncxantinuous  process,  and  because  the  traces  are  carried  ait  for  discrete 
values  of  the  parameters  used  to  define  the  situation  of  interest,  it  has 
been  of  prime  importance  to  develop  high  precisian  interpolation  and 
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extrapolation  techniques  to  allow  the  estimation  of  the  required  values  of 
the  variables  from  the  sets  of  discrete  parameter  values  available.  In 
addition  to  the  use  of  such  techniques  in  the  prediction  of  the  position 
and  development  parameter  values,  as  mentioned  in  the  preceding/  the 
techniques  are  critically  important  to  the  accurate  and  efficient 
determination  of  the  location  of  band  edges  in  ary  parameter  space 
(including  cutoff  values) ,  and  in  the  procedures  employed  in  the  automatic 
mapping  of  penumbral  structures. 


Uie  computer  program  can  operate  in  a  number  of  different  modes,  each 
designed  to  carry  out  one  of  the  well  defined  tasks  involved  in  the  various 
aspects  of  examining  and  mapping  the  penumbra.  Although  the  selection  and 
definition  of  the  individual  program  task  modes  is  under  interactive 
control,  the  execution  of  each  task  is  performed  completely  automatically 
by  the  ccnputer.  me  more  important  of  the  program  modes  of  operation  are 
described  in  the  following  sections. 


4.4.1  MAPPING:  me  mapping  facility  traces  the  edges  of  the  penumbral 
bands  by  means  of  a  continued  cycle  of  operations  in  which  two  parameters 
(zenith  and  azimuth,  or  latitude  and  longitude,  for  example)  are 
continuously  varied  while  the  matching  condition  is  continuously  applied. 
This  procedure  is  essentially  a  starting  point  calculation  in  the  two 
parameter  space,  in  which  the  locus  is  traced  of  paired  values  of 
parameters  for  which  grazing  occurs  at  the  lew  point  of  interest. 

The  domain  over  which  the  locus  (or  "sheet”,  as  Lemaitre  and  Vallarta 
termed  it)  contains  the  allowed-forbidden  transition  trajectories  is 
normally  limited  because  of  the  intruding  effects  of  other  trajectory 
minima.  Generally  some  lew  point  other  than  the  one  initially  considered, 
will  move  downward  with  continued  change  in  the  variable  parameters,  and  at 
sane  place  on  the  locus  this  second  low  point  will  lie  at  the  grazing 
height.  Beyond  this  the  newly  significant  low  point  will  cause  the 
trajectory  to  be  forbidden  so  that  now  the  further  "virtual"  extension  of 
the  locus  will  no  longer  define  the  allcwed-forbidden  transition  (although 
it  is  useful,  for  some  purposes,  to  map  such  virtual  penumbral  band  edges) . 
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The  transition  will  now,  for  some  finite  domain  within  the  two  parameter 
space,  be  defined  by  a  locus  associated  with  grazing  by  the  second  loop. 

The  whole  penumbral  structure  in  any  situation  is  constituted  of  sets  of 
such  intersecting  sheets. 

A  predictor-ccrrector  system  allows  the  program  to  "feel"  its  way 
along  any  given  sheet,  adjusting  the  size  of  its  steps  and  the  direction 
followed  according  to  its  predictive  success  in  the  preceding  steps.  In 
the  ocnputer  program  the  option  may  be  taken  of  tracing  the  system  of 
sheets,  automatically  transferring  fran  one  sheet  to  another  at  the  points 
of  intersection.  Alternatively,  the  entire  extent  of  a  sheet,  including 
its  virtual  extensions,  may  be  traced. 

Figure  4.6  shows  a  typical  example  of  the  kind  of  mapping  that  can  be 
carried  out.  This  example  shows  the  zenith-azimuth  position  of  the  various 
major  penumbral  structures  for  the  specified  conditions,  for  rigidity 
values  of  7.7,  14.0,  30.0  and  ®  GV.  Many  fine  bands  co-exist  here  too, 
lying  predominantly  close  to  and  parallel  to  those  bands  mapped.  The 
characteristics  of  the  major  structures  in  the  diagram  are  broadly  similar 
to  those  found  by  Lemaitre  and  Vallarta  (for  example,  see  figure  7  of 
Lsnaitre  and  Vallarta,  1936b) .  The  characteristic  distortion  of  the 
structures  in  the  north-east  is  produced  by  "folding"  associated  with  the 
loop  cone  effect  of  Cooke  and  Humble  (1970) .  The  7.7  GV  figure  displays 
some  very  significant  features,  in  particular  the  isolated  forbidden 
"island",  which  form  as  the  penumbral  bands  associated  with  the  higher 
rigidity  structures  separate  from  the  major  underlying  penumbral  edge  with 
diminishing  rigidity. 

Figure  4.7  shows  a  map  of  the  boundary  of  the  shadow  cone  for  30  GV 
particles,  as  traced  by  computer,  using  the  trajectory  parameterization 
technique  (refer  to  Chapter  1  for  definition  of  "shadow  cone") .  In  brief  - 
it  is  the  shadow  of  the  planet  (in  this  case  the  Earth)  for  particles  of 
the  given  rigidity.  It  is  interesting  to  note  that  the  shape  of  the  shadow 
is  not  a  siirple  circle  or  oval  as  might  intuitively  be  expected,  it  instead 
exhibits  a  cusp  at  its  highest  zenith  angle  excursion.  This  occurs  because 
of  the  development  of  an  intermediate  loop  (i.e.  between  the  "horizon"  and 
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FIGURE  4.6 

Penuntaral  structure  for  50,30,20,14,10  and  7.7  GV  at  the  point  in  the 
geomagnetic  field  with  geographic  latitude  10°  N  and  lcngituSe  270°,  and 
altitude  400  Ion.  The  magnetic  field  model  used  was  IGRF  updated  to  1980.0. 
The  letters  used  to  identify  portions  of  the  structures  far  which  Band 
Stability  Factors  have  been  calculated  (the  results  are  presented  in 
4.2).  The  reference  circles  lie  30°  apart  in  zenith,  and  the  radial  lines 
30°  in  azimuth. 
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the  arrival  site)  at  sane  particular  azimuth.  With  further  changing 
azimuth  this  loop  develops  to  the  point  where  it  intersects  the  surface  of 
the  Earth,  upon  which  it  then  defines  part  of  the  original  bounding  line. 


s 


FIGURE  4.7  Simple  shadow  cane 
calculated  for  30  GV  charged  cosmic  ray 
particles  approaching  a  site  with 
geographic  coordinates:  Latitude  10°, 
longitude  270°  (other  calculation 
parameters  are  the  same  as  used  for  the 
plots  in  figure  4.6). 


4.4.2  MIGRATION:  The  operation  of  migrating  a  penumbral  structure  is 
essentially  similar  to  mapping.  In  this  case  the  objective  is  to  track  the 
position  of  a  given  penumbral  band  edge  as  it  moves  in  response  to  change 
in  any  of  the  input  parameters,  rather  than  to  acnpletely  ascertain  the 
form  of  the  penumbral  structure  in  the  domain  of  interest.  Changes  in  band 
edge  position  occurring  in  response  to  changes  in  a  number  of  parameters 
can  be  established  by  migrating  the  band  edge  either  one  parameter  at  a 
time  or  alternatively  by  simultaneous  stepped  proportional  change  in  all 
parameters. 


Table  4.1  shews  an  example  of  multiple  parameter  migration.  In  this 
example,  the  band  edge  which  exists  for  the  conditions  shewn  in  step  1  is 
to  be  traced,  and  the  longitude  at  which  the  band  will  lie  at  for  the  new 
situation  listed  in  step  7,  determined.  During  this  migration  each  of  the 
parameters  other  than  longitude  is  being  changed  proportionally.  longitude 
value  is,  in  this  case,  being  allcwed  to  "float"  -  an  appropriate  value 
being  determined  for  each  of  the  steps. 

In  single  parameter  migration,  one  parameter  is  being  changed  while 
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the  other,  nominated,  parameter  is  allowed  to  float.  In  practice,  for  the 
case  of  single  parameter  migration,  it  has  proved  convenient  to  use  zenith 
angle  as  the  dependent  (or  floating)  parameter,  so  that  the  position  of  the 
band  edge  is  followed  in  zenith  angle  as  each  parameter  is  varied. 


TABLE  4.1 

Example  of  multiple  parameter  "migration".  The  penumbral  band  edge 
which  lies  at  the  location  given  in  step  1  is  traced  to  find  at  what 
longitude  it  would  lie  for  the  parameter  values  listed  in  step  7  (note  that 
the  longitude  value,  initially  defined,  is  allowed  to  "float",  in  order  to 
find  what  new  value  of  longitude  would  plaoe  the  same  band  edge  in  the  new 
location  and  direction.  Angles  are  expressed  in  degrees. 


STEP 

LATITUDE 

LONGITUDE 

ZENITH 

AZIMUTH 

RIGIDITY 

ALT1TULE 

1 

7.5 

270.0 

0.0 

0.0 

10.0 

400.0 

2 

7.7 

(268.5) 

0.4 

24.1 

10.1 

395.6 

3 

8.8 

(266.1) 

1.4 

76.3 

10.3 

386.1 

4 

8.7 

(262.4) 

2.4 

131.3 

10.5 

376.1 

5 

9.2 

(255.8) 

3.4 

186.3 

10.7 

366.1 

6 

9.7 

(247.3) 

4.4 

241.3 

10.9 

356.1 

7 

10.0 

(243.1) 

5.0 

275.0 

11.0 

350.0 

4.4.3  PENUMBRAL  BAND  STABILITY:  The  stability  of  ary  given  portion  of  a 
penumbral  band  edge  can  very  siirply  be  expressed  in  terms  of  the  "Band 
Shift  Factor",  BSF  -  SR/ SR,  and  determined  by  computer  as  BSF  •  (6R/6Ri) 
t Pi  o  where  SR  is  the  change  in  rigidity  R  required  to  maintain  grazing  in 
the  low  point  responsible  for  the  penumbral  feature,  for  a  change  JPi  in 
the  1th  parameter  P^.  The  BSF,  being  the  measure  of  the  rate  at  which  the 
band  edge  shifts  with  variation  in  the  independent  parameter,  is 
essentially  a  factor  describing  the  stability  of  the  structure.  The 
independent  parameters  used  to  establish  the  BSFs  could  be  extended  to 
include  quantitative  measures  of  individual  internal  or  external  field 
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sources,  internal  calculation  parameters,  or  any  other  "real"  or  artificial 
parameter  involved  in  the  calculations.  The  ccmputer  program  has  been  set 
up  to  produce,  on  demand,  a  full  set  of  BSFs  pertaining  to  a  nominated  band 
edge,  for  the  basic  set  of  input  parameters. 


TABLE  4.2 

Values  of  Band  Shift  Factor,  a  measure  of  penumbral  band  edge  stability, 
for  the  structures  indicated  alphabetically  in  figure  4.6. 


I 

IBand  shift  factor  value  for  different  structures  at  Identified  locations 


Independent 

Parameter 

1  Absolute  1 
j  Value  of  i 

I  Parameter! 

A  B  C  0  E  F 

G  H 

Unit 

Latitude 

T  1 

1  10°  1 

2.6772  0.3131  O.UOIB  0.0319  0.0107  0.0044 

0.4660  -0.0079 

GV/* 

Longitude 

|  270°  1 

-0.1510  0.1993  -0.0010  -0.0600  -0.0510  -0.0529 

-0.0507  -0.0002 

r,v/* 

Geocentric  alt. 

1  1 
|  400  km  | 

-0.0395  -0.00%  -0.0033  -0.0009  -0.0029  -0.0031 

0.0095  -0.0523 

GV/kei 

Grating  height 

1  1 

1  30  km  I 

0.0247  0.0069  0.0004  0.0033  -0.0001  0.0001 

-0.0542  -0.0001 

GV/km 

Zenith  angle 

1  1 

I  various*  1 

-0.4212  -0.3299  -0.2235  -0.2827  -0.2119  -0.2290 

0.0038  -0.0814 

GV/* 

Azimuth  angle 

i  i 
I  various*  1 
l*see  flg.ll 

-0.2932  -0.0006  -0.0663  -0.1231  -0.0837  -0.0816 

-0.0128  -0.0263 

GV/* 

The  dependence  of  the  change  in  penumbral  band  edge  position  in  any 
parameter  Pj  for  change  in  any  other  parameter  Pj  can  be  obtained  by 
suitable  combination  of  the  BSF's,  viz. 


3P2  3P2  3R 


Similarly,  the  stability  of  the  entire  penumbral  band  can  be  expressed 
by  suitable  combination  of  BSFs  to  form  the  Band  Growth  Factor  (BGF) : 
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BCF 


(3R„/3P  -  3R  /  3P ) 

_ t _ ± _ _  x  100 

RU  ~  RL 


X 


per  unit  change  in  the  independent  parameter  P,  where  Py  and  Rl  are  the 
rigidities  defining  the  upper  and  lower  edges  of  the  penumbral  hand.  The 
BGF  is  a  measure  of  the  percentage  rate  of  change  in  width  of  a  penumbral 
band  with  change  in  any  given  parameter  for  a  given  set  of  conditions. 

Table  4.2  presents  calculated  values  of  the  Band  Shift  Factors 
pertaining  to  eight  points  an  the  penumbral  structures  presented  in  figure 
4.6,  where  the  points  are  labeled  alphabetically.  Note  the  tendency  far 
the  stability  to  be  greater  for  larger  scale,  short  range,  structures,  and 
fear  higher  rigidities  (a  high  BSF  implies  high  stability) . 


4.5  OONCUJSICNS 

The  trajectory  parameterization  technique  offers  a  very  powerful  means 
of  determining  the  characteristics  of  the  penumbra  in  any  particular 
situation.  Use  of  the  technique  has  revealed  seme  very  interesting,  and 
seme  unexpected,  information  about  the  penumbra.  For  example,  the 
penumbral  "islands"  displayed  in  figure  4.6  were  not  known  to  exist 
previously.  Use  of  the  migration  and  stability  techniques  show  these 
islands  to  have  interesting  properties,  in  particular,  a  great  sensitivity 
to  rigidity  and  field  model  (for  example,  the  7.7  GV  penumbral  island, 
centered  at  76°  zenith  angle  and  279°  azimuth  in  figure  4.6,  disappears 
when  the  rigidity  is  reduced  by  0.085  GV,  and  does  not  exist  at  7.7  GV  when 
the  mapping  is  carried  out  using  the  IGRF  1975  geomagnetic  field) . 

In  the  following  chapters  the  parameterization  technique  is  used  to 
investigate  other  aspects  of  the  penumbra  within  the  geomagnetic  field. 
These  techniques  would  apply  equally  well  to  the  magnetic  field  of  Jupiter, 
for  example,  but  to  date  have  not  yet  been  applied  to  fields  other  than 
that  of  the  Earth. 
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CHAPTER  5:  TOE  HIGH  ZENITH  ANGLE  LIMITS  OF  COSMIC  RAY  ACCESS  TO  AN  EARTH 
SATELLITE 

5.1  INTOODUCTICN 

A  knowledge  Is  often  required  of  the  highest  zenith  angle  far  which  a 
charged  particle  of  given  rigidity  has  access  to  an  earth  satellite  lying 
in  a  particular  orbital  configuration. 

Because  of  the  complexities  of  the  trajectories  by  which  below-line-of 
sight  access  to  a  satellite  occurs  it  has  proven  difficult  to  determine 
systematically  whether  absolute  limits  exist  beyond  which  access  is 
forbidden.  The  general  approach  to  this  problem  has  been  to  directly  test 
for  access  by  carrying  out  a  large  volume  of  trajectory  traces.  This 
method,  however,  is  not  completely  satisfactory  because  the  absence  of 
detected  access  beyond  a  given  angle  does  not  constitute  a  satisfactory 
proof  that  access  could  not  occur  (either  via  undetected  allowed 
trajectories  that  exist  under  the  geomagnetic  field  conditions  simulated, 
or  under  even  slight  perturbations  away  from  these  conditions) . 

A  technique  has  been  developed,  based  on  the  trajectory 
parameterization  technique,  which  allows  a  very  much  more  definitive 
insight  into  the  problem  of  establishing  the  high  zenith  angle  limits. 


5.2  DISCUSSION 

The  new  method  involves  the  examination  of  the  form  of  trajectories 
along  which  particles  approach  a  satellite,  and  consideration  of  whether 
these  orbits  can  or  cannot  provide  access  for  charged  particles  entering 
from  outside  the  geomagnetic  field.  (In  line  with  ordinary  trajectory 
tracing  practice  the  traces  are  carried  out  in  the  reverse,  negative-time, 
direction  -  out  away  from  the  satellite.) 

Before  discussing  the  technique  further  it  is  appropriate  to  briefly 
review  the  question  of  access  to  a  point  in  an  axially  symmetric 
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representation  of  the  geccagnetic  field  (more  fully  rij  scussed  in  Chapter 

1) .  In  a  simple,  axially  symmetric  field,  three  types  of  solid  angle 
region  may  be  distinguished  relative  to  any  point  in  the  field  -  the  "main 
allowed"  cone,  the  "Stormer"  cone,  and  the  penumbra. 

For  the  purpose  of  the  work  reported  here  it  has  proven  practical  and 
useful  to  distinguish  four  different  kinds  of  solid  angle  regions  relative 
to  a  given  point  in  the  "real"  field.  These  are  the  regions  within  which: 

1)  access  is  forbidden  because  of  short  range  intersection  with  the 
Earth  (as  detected  in  the  outward  trajectory  trace) .  These  trajectories 
correspond  to  the  short  range  "shadow"  orbits  of  Lemaitre  and  Vallarta. 

2)  access  from  outside  the  geomagnetic  field  either  does,  or 
possibly  ocaild,  occur  along  aperiodic  or  unstable  quas i-per iodic  orbits. 

3)  entry  would  only  be  by  aperiodic  of  quas  i-per  iodic  orbits  that 
intersect  the  surface  of  the  earth. 

4)  approach  oculd  occur  via  quasi-periodic  orbits  that  have  no 
detected  intersection  with  the  atmosphere.  Such  trajectories  can  extend 
for  very  large  distances  in  the  field  without  any  detected  field  escape  or 
earth  intersection.  It  is  impractical  (or  even  sometimes  impossible)  to 
trace  these  to  some  definite  end,  so  the  traces  are  terminated  at  a 
nominated  path  length.  Thus  the  matter  of  access  via  such  trajectories 
remains  unresolved  (the  question  of  the  precision  of  very  long  trajectory 
traces  necessarily  enters  here,  but  will  not  be  discussed) . 

Figures  5.1  through  5.4  illustrate  trajectories  associated  with  the 
four  kinds  of  region.  In  the  trajectory  parameterization  technique  used  in 
this  study  it  is  important  to  continue  the  trajectory  traoes  on  beyond  the 
first  intersection  with  the  surface  of  the  Earth.  This  is  an  artifact 
necessary  to  the  operation  of  the  technique  (and  in  the  identification  of 
trajectory  form) ,  and  is  not  held  to  represent  any  kind  of  reality.  Thus, 
in  figures  5.1  and  5.3,  the  trajectories  are  seen  to  pass  through  the 
surface  of  the  solid  Earth  (shown  by  a  dashed  line) . 
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FIGURE  5.1 


Trajectory  (type  (1))  displaying  short  range  "shadcw4'  intersection 
with  the  atmosphere.  The  trajectory,  which  was  continued  on  past  the 
intersection,  indicated  that  escape  from  the  field  (or  entry  from  outside 
the  field)  would  otherwise  have  been  possible. 

[Trajectory  arrival  point  parameters:  latitude  ■  17.5°  N,  longitude  «  260° 
E,  zenith  -  145°,  aziruth  -  270°,  rigidity  -  7.5  GV,  Altitude  -  400  to, 
grazing  height  *  30  to;  geomagnetic  field  IGRF80] . 


o  -0.01 


Satoliite 


FIGURE  5.2 


LONGITUDE 


Trajectory  of  type  (2)  -  along  vhich  the  point  in  the  field  is 
accessible  via  aperiodic  orbit. 

[Trajectory  parameters:  latitude  ■  10°,  longitude  •  270°,  zenith  »  144.5°, 
azimuth  •  278°,  rigidity  -  7.7,  altitude  -  400  tan,  grazing  height  -  30  km; 
geomagnetic  field  IGRF80] . 
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FIGURE  5.3 
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360/0* 


Type  (3)  trajectory  -  shewing  regular  quasi-bound  periodic  farm, 
intersecting  the  atmosphere. 

[Trajectory  parameters:  latitude  -  7.5°  N,  longitude  *  45°,  zenith  •  145°, 
azimuth  ■  270°,  rigidity  ■  7.5  GV,  altitude  ■  400  km,  grazing  height  *  so 
km;  geomagnetic  field  IGRF80] . 
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Preliminary  investigations  have  been  carried  out  to  elucidate  the  form 
and  properties  of  the  regions  containing  the  trajectories  of  the  four 
kinds.  As  a  first  step  in  this  study,  the  regions  containing  trajectories 
of  type  (4)  are  mapped,  either  in  z enith-az imuth ,  or  in  latitude-longitude, 
spaoe.  An  example  of  zenith-azimuth  mapping,  as  given  by  Cooke,  1982,  is 
reproduced  in  part  as  figure  5.5.  She  high  zenith  angle  limits  are  set  by 
the  structure  visible  in  the  west  at  *145°  zenith  angle.  This  boundary  is 
formed  by  short  range  shadow  orbits. 

When  mapped  in  latitude  and  longitude  the  regions  typically  have  the 
form  shewn  in  figure  5.6.  At  the  given  zenith  and  azimuth  (150°  and  270° 
respectively)  access  by  positively  charged  7  GeV  particles  to  a  satellite 
in  a  400  km  geooentric  orbit  is  explicitly  forbidden  at  all  points  on  the 
Earth's  surface  shown  by  the  hatching.  The  boundary  of  the  central  region 
(within  which  access  is  not  explicitly  forbidden  as  a  result  of  1st  and  2nd 
loop  Earth  intersections)  is  in  part  formed  by  the  edge  of  the  region 
associated  with  1st  loop  intersections,  and  in  part  by  the  edge  of  the 
region  associated  with  the  2nd  loop  intersections.  As  rigidity  varies  the 
region  defined  by  the  first  order  shadow  orbits  changes  shape  as  shown  in 
figure  5.7.  The  2nd  order  shadow  limits  are  more  stable,  and  change  very 
little  over  this  rigidity  range.  A  charge  in  zenith  typically  produces  the 
shift  shewn  in  Figure  5.8,  and  change  in  azimuth  the  shift  to  the  positions 
shewn  in  figure  5.9  and  5.10. 

The  Band  Shift  Factor  (see  Chapter  4)  can  be  used  to  quantitatively 
express  the  stability  of  these  structures  (see  Table  5.1) .  It  is 
interesting  to  note  the  relatively  great  stability  of  the  central  second 
order  structure,  which  lies  along  the  geomagnetic  equator. 

By  examination  of  the  form  and  position  of  the  shadow  orbit  defined 
regions  it  is  possible  to  discount  the  possibility  of  high  zenith  angle 
entry  at  azimuths  well  away  from  the  west,  and  to  rigidities  less  than  a 
given  value  (which  value  depends  an  the  zenith  angle  considered,  e.g.  entry 
is  forbidden  for  rigidity  >  9  GV  at  zenith  angle  of  150°) . 
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FIGURE  5.5 

Position  of  allowed  and  forbidden  structures  mapped  in  zenith  and 
azimuth  space  for  the  indicated  rigidity  values,  at  the  point  10°  N 
latitude,  270  E  longitude  (IGRF80  field) .  The  shadow  structure  limiting 
the  high  zenith  angle  entry  in  the  west,  for  the  7.7  GV  rigidity  value,  is 
indicated  by  an  asterisk. 


FICKJRE  5*6 

latitude-longitude  map  showing  zones  of  aooess  under  the  indicated 
conditions.  The  Band  Shift  Factors  have  been  calculated  for  points  A,  B, 
and  C,  and  are  presented  in  Table  5.1.  The  shaded  region  is  inaccessible 
due  to  intersection  of  1st  and/or  2nd  loops  with  the  atmosphere. 


Latitude  vs  Longitude 

Zenith  *  146.000  Azimuth  «  270.000  Rigidity  --- 

Altitude  =  100.000  Reentry  =  30.000  Year  =  ; 
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1st  LOOP-FORBIDDDEN 


2nd  LOOP-FORBIDDEN 
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FIGURE  5.8 


Latitude-longitude  map  for  the  conditions  noted,  showing  categorized 
zones  of  aooess  (this  figure  corresponds  to  the  sane  set  of  conditions 
under  vtaich  figure  5.6  was  produced,  except  that  it  pertains  to  145°  zenith 
rather  than  the  150°  zenith  of  Figure  5.6) . 


latitude 


latitude  vs  Longitude 

Zenith  =  150.000  Azimuth  =  315.000  Rigidity  =  7.000 

Altitude  =  400.000  Reentry  =  30.000  Year  =  1980.000 
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2nd  LOOP-FORBIDDEN 
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FIGURE  5.9 


Latitude- longitude  map  drawn  for  same  conditions  as  for  Figure  6 
except  that  the  azimuth  is  new  225°,  rather  than  270°  east  of  north. 


Latitude  vs  Longitude 

Zenith  -  150.000  AzunuLh  -  225.000  Rigidity  -•<  7  000 

Altitude  -  400.000  Reentry  =  30.000  Year  *  1080.000 
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2nd  LOOP-FORBIDDEN 


LONGITUDE 


FIGURE  5.10 


5.6, 


Latitude-longitude  nap  drawn  for  the 
except  that  the  azimuth  is  now  315°, 


sane  conditions  as  for  figure 
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TABLE  5.1 


Band  Stability  Factor  S'R/S'P,  where  R  is  rigidity,  and  P  is  the 
parameter  in  question,  for  the  points  A,  B,  and  C  in  Figure  5.6.  Large  BSF 
values  are  associated  with  relatively  stable  structures. 


POINT 

LAT 

LON 

ZE 

AZ 

ALT 

GRAZ 

A 

-0.595 

-0.043 

-0.174 

-0.054 

0.015 

-0.018 

B 

-0.750 

-0.030 

0.236 

-0.247 

-0.014 

0.013 

C 

0.109 

-0.004 

-0.245 

0.032 

0.019 

-0.023 

Within  the  central  region  exist  fine  allcwed-forbidden  structures 
associated  with  3rd,  4th,  and  so  on,  loop  intersections.  The  mapping  of 
these  structures  has  been  found  impractical  because  of  the  large  amount  of 
computer  time  required.  In  any  case  the  detailed  form  of  the  structures 
depend  significantly  on  the  choice  of  geomagnetic  field  model  used  in  the 
calculations  (whereas  the  structures  associated  with  1st  and  2nd  loop 
intersections  are  relatively  stable) .  It  is  therefore  necessary  to  use 
some  other  means  of  determining  what  access  is  possible  into  this  region. 

Examination  of  the  form  of  the  trajectories  associated  with  approach 
to  points  within  the  central  region  allows  the  question  of  access  from 
points  outside  the  geomagnetic  field  to  be  pursued,  and  permits  general 
conclusions  to  be  drawn  about  whether  access  is  possible  or  impossible  in 
given  situations.  This  examination  has  been  done  by  means  of  sets  of 
trajectory  traces  applying  to  directions  spaced  in  latitude  and  longitude 
(typically  at  5°  intervals  in  latitude,  and  30°  in  longitude) .  The 
trajectories  are  displayed  on  a  video  terminal  and  visually  categorized 
into  types  2,3,  or  4.  Generally  type  3  and  4  orbits  are  readily 
distinguished,  although  there  is  a  less  distinct  difference  between  type  3 
and  type  2  orbits.  The  question  arises  of  when  short  range  quasi- 
periodicity  turns  into  a  regular  periodicity.  Hot  withstanding  this  point, 
meaningful  boundaries  can  be  established.  Figure  5.8  shews  the  location  of 
the  boundaries  for  this  set  of  conditions. 


In  general  it  is  found  that  trajectories  originating  close  to  the 
geomagnetic  equator  shew  the  most  simple  and  regular  periodicity  (see 
figures  5.3  and  5.4} .  Further  away  to  either  side  the  trajectories  beocne 
more  complex  in  their  periodic  form,  and  further  away  again  degenerate  into 
unstable  quas i-per iodic  or  aperiodic  forms  (see  figure  5.2),  along  which 
access  may  occur  from  outside  the  field  (the  short  range  shadow  structures 
tend  to  be  responsible  for  preventing  access  along  what  otherwise  may  be 
allowed  trajectories) .  It  is  possible  to  anticipate  whether  trajectories 
of  type  2  occur  before  the  shadow  edges  close  off  access,  and  it  is 
therefore  found  unnecessary  to  actually  locate  an  allowed  trajectory  in 
order  to  anticipate  possible  access  via  trajectories  of  this  type. 

The  offset  of  the  earth's  dipole  gives  rise  to  a  strong  longitudinal 
asymmetry  in  the  form  of  the  regions  (see  figure  5.7,  for  example) .  There 
is  a  strong  tendency  for  type  4  trajectories  to  be  normally  associated  with 
aooess  to  points  at  longitudes  close  to  300°,  whereas  type  3  predominate 
over  the  rest  of  the  longitude  range.  Because  the  1st  order  shadow 
patterns  jure  mainly  (for  the  conditions  of  interest)  centered  within  the 
range  90-180°  longitude,  these  regions  can  be  regarded  as  being  populated 
largely  by  type  3  trajectories,  with  type  2  orbits  (if  they  occur  at  all  in 
a  given  situation)  lying  at  points  removed  from  the  geomagnetic  equator. 

In  general,  as  rigidity  decreases,  the  size  of  the  mapped  regions 
increase.  At  the  same  time,  however,  the  trajectory  configurations  evolve 
into  a  much  more  tightly  bound  periodic  form,  of  types  3  and  4.  Particles 
could  conceivably  populate  such  trajectories  -  possibly  primary  cosmic  rays 
under  time  varying  conditions,  or  splash  albedo  particles,  for  exairple. 

5.3  OONOUSICN 

In  extending  this  preliminary  investigation  it  will  be  necessary  to 
mate  an  examination  of  the  zenith,  azimuth,  rigidity,  and  altitude 
dependence  of  the  mapped  regions.  In  this  way  a  more  ccnpl ete  picture 
could  be  drawn  up  of  the  zenith  angle  limits  of  aooess  -  limits  to  the 
regions  containing  trajectories  associated  with  aooess  of  the  kinds 
possible  and  impossible,  probable  and  improbable. 
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CHAPTER  6:  THE  MAIN  CONE  IN  THE  REAL  MAGNETIC  FIELD 

6.1  3NTRODUCTICN 

The  "main  ccrne"  is  the  boundary  of  the  fully  unconditionally  allowed 
cone  within  which  oosmic  rays  of  a  given  rigidity  approach  a  point  within  a 
planetary  magnetic  field.  It  is  constituted  in  part  by  trajectories  which 
are  asymptotic  to  the  simplest  bound  periodic  orbits,  and  in  part  by 
trajectories  which  graze  the  surface  of  the  planet.  In  this  chapter  the 
phenomenology  of  the  former  of  these  two  types  of  trajectories  is 
discussed,  and  a  technique  is  described  by  which  main  cone  cutoffs  may  be 
calculated  by  computer  to  very  high  precisian. 

Lemaitre  and  Vallarta  investigated  the  properties  of  the  main  cone  as 
it  applies  to  a  simple  dipole  representation  of  the  Earth's  magnetic  field. 
The  currently  reported  research  has  shewn  that  a  corresponding  true  main 
cone  also  exists  at  locations  in  the  real  geomagnetic  field.  The  form  of 
the  main  cone  in  the  real  field  is  very  similar  to  that  found  by  Lemaitre 
and  Vallarta  to  occur  in  the  simple  dipole  field  representation. 

The  facility  of  determining  the  main  cane  structure  in  the  real  field 
is  valuable,  because  often  a  knowledge  is  required  of  the  rigidity  limit 
above  which  unconditional  oosmic  ray  entry  could  be  expected  to  exist  in 
specific  situations  under  real  field  conditions.  The  normal  unmodified 
iterative  trajectory  trace  method  of  penunibral  mapping,  relying  as  it  does 
on  the  detection  of  forbidden  trajectories,  cannot  yield  a  reliable 
estimate  of  this  upper  rigidity  limit. 

6.2  DISCUSSION 

Figure  6.1  shows  a  typical  main  cane  trajectory,  which  displays 
characteristic  regular  equator  crossings.  Note  that  the  trajectory,  like 
others  of  its  kind,  lies  essentially  within  a  constant  altitude  shell, 
disposed  at  an  essentially  constant  distance  fretn  the  Earth's  equivalent 
dipole.  The  dependence  of  the  position  of  the  shell,  as  determined  from 
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FIGURE  6.1 

Typical  "read."  field  asynptotic  quasi-bound  periodic  orbit,  which  lies 
very  close  to  the  main  cone  edge  (trajectory  starting  point  parameters: 
geographic  latitude  10°  N,  longitude  270°  E,  zenith  angle  7.9368°,  azimuth 
270°,  rigidity  10  GV,  400  km  geocentric  altitude?  IGRF80  field 
representation) .  The  position  of  the  trajectory  is  plotted  as  a  function 
of  geographic  latitude  and  longitude,  and  geooentric  altitude. 


trajectories  calculated  to  exist  in  a  real  field  representation,  an 
rigidity  corresponds  closely  to  the  Storxner  radius  for  the  given  rigidity 
(the  radius  at  which  a  charges  particle  traces  a  circular  orbit  about  the 
Earth's  equivalent  dipole).  Figure  6.2  shews  the  dependence  of  the 
position  of  the  shell  as  a  function  of  rigidity,  both  as  found  to  exist  in 
a  real  field  representation,  and  as  predicted  by  the  Stonier  relation 

radius  *=  J (59.6/rigidity)  for  the  Earth's  field 

It  is  useful  to  appreciate  that  main  cons  cutoff  values  of  more  than  about 
60  GV  cannot  exist  in  the  real  field  because  the.  shall  containing  the  bound 
periodic  orbit  would  have  to  lie  within  the  solid  Earth.  Neither  could 
main  cane  cutoff  values  below  about  0.5  GV  bo  found  in  the  real  field 
because,  in  this  case,  the  shell  containing  the  bound  periodic  orbit  would 
have  to  lie  at  a  radius  of  more  than  10  Earth  radii  from  the  Earth's 
equivalent  dipole  -  which  would  place  the  shell  outside  the  magnetospheric 
boundary  an  the  "unwind"  (i.e.  solar  wind)  side  of  the  field. 

Trajectories,  for  specific  rigidity  values,  traced  in  the  "negative 
time"  direction,  from  the  point  of  interest  in  the  field  cutward  through 
the  magnetosphere,  cut  through  the  shells  corresponding  to  lower  rigidity 
values  with  a  positive  slope.  In  particular,  at  the  point  along  the 
trajectory  where  the  magnetic  equator  is  crossed,  the  slope  of  the 
trajectory  at  those  points  is  positive.  On  the  other  hand,  it  is  found 
that  the  slopes  at  the  corresponding  equator  crossings  associated  with 
rigidity  values  lower  than  the  main  com  cutoff  are  normally  negative. 

The  present  method  of  determining  precise  main  cone  cutoff  values 
utilizes  this  fact,  and  the  corputer  program  which  evaluates  the  cutoffs 
locates  the  equator  crossings  and  determines  their  slopes.  Because, 
however,  trajectories  pertaining  to  rigidities  well  removed  (mainly  higher) 
from  the  main  cone  cutoff  may  have  few  of  no  equator  crossings,  it  is 
necessary  first  to  derive  an  initial  estimate  of  the  main  cone  cutoff 
before  refining  it  by  use  of  the  equator  crossing  slopes. 

To  this  end  the  Stormer  equation  is  invoked .  The  resulting  rigidity 
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trajectories  calculated  to  exist  in  a  real  field  representation,  on 
rigidity  (corresponds  closely  to  the  Stormer  radius  for  the  given  rigidity 
(the  radius  at  which  a  charges  particle  traces  a  circular  orbit  about  the 
Earth's  equivalent  dipole).  Figure  6.2  shows  the  dependence  of  the 
position  of  the  shell  as  a  function  of  rigidity,  both  as  found  to  exist  in 
a  real  field  representation,  and  as  predicted  by  the  Stormer  relation 

radius  *=  J (59.6/rigidity)  for  the  Earth's  field 

It  is  useful  to  appreciate  that  main  cone  cutoff  values  of  more  than  about 
60  GV  cannot  exist  in  the  real  field  because  the  shell  containing  the  bound 
periodic  orbit  would  have  to  lie  within  the  solid  Earth.  Neither  could 
main  cone  cutoff  values  below  about  0.5  GV  be  found  in  the  real  field 
because,  in  this  case,  the  shell  containing  the  bound  periodic  orbit  would 
have  to  lie  at  a  radius  of  more  than  10  Earth  radii  from  the  Earth's 
equivalent  dipole  -  which  would  place  the  shell  outside  the  magnetospheric 
boundary  on  the  "upwind"  (i.e.  solar  wind)  side  of  the  field. 

Trajectories,  for  specific  rigidity  values,  traced  in  the  "negative 
time"  direction,  from  the  point  of  interest  in  the  field  outward  through 
the  magnetosphere,  cut  through  the  shells  corresponding  to  lower  rigidity 
values  with  a  positive  slope.  In  particular,  at  the  point  along  the 
trajectory  where  the  magnetic  equator  is  crossed,  the  slope  of  the 
trajectory  at  those  points  is  positive.  On  the  other  hand,  it  is  found 
that  the  slopes  at  the  corresponding  equator  crossings  associated  with 
rigidity  values  lower  than  the  main  cone  cutoff  sue  normally  negative. 

The  present  method  of  determining  precise  main  cone  cutoff  values 
utilizes  this  fact,  and  the  computer  program  which  evaluates  the  cutoffs 
locates  the  equator  crossings  and  determines  their  slopes.  Because, 
however,  trajectories  pertaining  to  rigidities  well  removed  (mainly  higher) 
from  the  main  cone  cutoff  may  have  few  of  no  equator  crossings,  it  is 
necessary  first  to  derive  an  initial  estimate  of  the  main  cone  cutoff 
before  refining  it  by  use  of  the  equator  crossing  slopes. 

To  this  end  the  Stormer  equation  is  invoked.  The  resulting  rigidity 
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value  is  a  very  rough  estimate  -  to  within  only  about  40%.  This  estimate 
is  still  not  good  enough  to  start  the  cutoff  refining  procedures,  so  in 
order  to  inprove  the  estimate,  two  trajectories  are  traced  which,  although 
only  a  small  incremental  rigidity  value  apart,  are  clearly  above  the  actual 
main  cone  cutoff  value  (the  chosen  rigidity  values  are,  in  fact,  related  to 
the  Stormer  estimate  multiplied  by  1.5) .  The  values  of  the  slope  of  these 
trajectories  as  they  pass  through  the  shell  altitude  at  which  a  bound 
periodic  orbit  could  exist  at  their  particular  rigidity  values  are  then 
used  to  produce  a  more  precise  estimate  of  the  main  cone  cutoff  value. 

These  estimates  are  accurate  to  within  about  5%  -  quite  accurate  enough  to 
start  the  cutoff  value  refining  process. 

A  pair  of  trajectories  are  traced  for  this  rigidity  value  plus  and 
minus  a  small  rigidity  increment.  This  pair  of  trajectories  are  found 
generally  to  contain  three  or  more  equator  crossings.  The  values  of  the 
slopes  of  the  third  equator  crossings  are  noted,  and  the  dependence  of  the 
slope  on  rigidity  is  deduced.  The  value  of  the  rigidity  for  which  the 
slope  would  be  zero  is  then  ccrputed. 

A  new  incremental  pair  of  trajectories  are  traced  for  this  new 
rigidity  value  (the  increment  is  decreased  by  approximately  an  order  of 
magnitude  each  time  a  new  incremental  pair  of  trajectories  is  used)  and  the 
slope  at  the  fourth  equator  crossing  vised  to  produce  a  new,  refined, 
estimate  of  the  main  cone  cutoff.  Successive  incremental  pairs  of 
trajectories  are  traced  and  the  position  of  successive  equator  crossings 
used  to  progressively  adjust  the  rigidity  value. 

This  series  of  rigidity  estimates  asymptotically  approaches  the  "true" 
main  cone  value.  It  is  possible,  in  most  situations,  to  trace  main  cane 
edge  trajectories  for  many  equator  crossings  after  the  refining  process  has 
operated  for  a  few  cycles  of  operation  -  ten  crossings  is  generally 
possible.  The  precision  achieved  in  the  main  cone  cutoff  estimate  is 
correspondingly  great,  if  in  the  limit,  unreal,  due  to  limitations  on  the 
precision  of  modeling  the  geomagnetic  field.  For  most  purposes  it  has  been 
found  sufficient  to  aim  for  a  precision  of  0.001  GV,  which  involves  tracing 
between  ten  and  twenty  trajectories  for  between  5  art!  7  equator  crossings 


FIGURE  6.3 

Contours  of  mail  cutoff  rigidity  as  a  function  of  geographic  latitude 
and  longitude  at  1  GV  intervals  in  both  hemispheres. 
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FIGURE  6.5 

As  for  figure  6.4,  for  7.7  GV. 
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case  in  penumbral  trajectories,  but  only  to  the  particular  equator  crossing 
targeted  at  the  given  stage  in  the  ocnputaticn  procedure. 

Table  6.1  presents  a  list  of  rigidity  values  which  were  produced  in  an 
actual  main  cutoff  determination ,  which  shows  how  the  values  very  quickly 
asynptote  to  the  final  high  precision  value. 

TABLE  6.1 

A  typical  list  of  rigidity  values  associated  with  a  main  ocne  cutoff 
determination,  showing  the  progressive  refinement  of  the  estimates,  as 
successive  equator  crossings  are  examined.  In  this  case  the  precision 
required  of  the  final  cutoff  value  was  0.001  GV. 

EQUATOR  CROSSING 

2 

3 

4 

5 

6 

7 


It  should  be  noted  that  the  procedure  is  unaffected  by  the  form  of  the 
trajectory  in  the  area  of  the  starting  point.  Loops  in  the  magnetic  field 
line  bundle  in  this  region,  via  which  entry  into  the  asymptotic  trajectory 
must  be  made,  are  responsible,  as  discussed  by  Cooke  and  Huiribla  (1970) ,  far 
the  detailed  "folded"  structure  in  the  form  of  the  main  cone. 

Nevertheless,  the  sain  cone  edge  may  be  located  regardless  of  the 
cccplexity  of  this  part  of  the  orbit  and  for  any  arrival  direction  at  the 
site  of  interest. 


RIGIDITY 

14.527  (initial  Stormer-based  estimate) 
10.424  (refined  1st  estimate) 

9.924243 

9.931371 

9.927419 

9.927419 

9.925816 

9.925418 

9.925425  (final  main  cutoff  value) 
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The  trajectory  parameterization  method  of  deriving  the  main  cone 
cutoff  values  is  completely  automatic,  requiring  only  initial  insertion  of 
parameters  specifying  the  location  and  direction  of  interest  in  the 
magnetic  field.  Examples  of  the  results  obtainable  by  the  method  are  given 
in  figures  6.4  and  6.5. 

The  results  of  this  method  agree  with  the  main  cutoffs  derived  by  Shea 
and  Smart  using  a  procedure  that  enables  them  to  detect  the  main  cutoff 
while  the  standard  iterative  penumbral  mapping  is  proceeding.  Their  method 
essentially  relies  an  the  same  phenomenology  as  the  parameterization 
method. 
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CHAPTER  7:  CONCLUSIONS 


7.1  USE  OF  THE  "TRAJECTORY  PARAMETERIZATION"  PENUMBRAL  MAPPING  TECHNIQUE 

A  range  of  techniques  are  new  available  to  allcw  the  detailed  napping 
of  cosmic  ray  access  to  locations  within  planetary  nagnetospheres. 

Oanputer  programs  embodying  these  techniques  have  been  developed  that  traoe 
the  main  cutoff  value  in  any  direction,  traoe  penumbral  band  structure  over 
zenith-azimuth  space,  or  longitude-latitude  space,  and  establish  the 
stability  of  these  in  any  direction.  The  assembly  of  these  directional 
cutoffs  into  integrated  cutoff  distributed  functions  is  also  now  a  routine 
operation.  The  development  of  a  further  technique  was  oentnenoed  but  not 
completed  -  a  system  by  which  very  finely  detailed  penumbral  structure  in 
any  direction  may  be  automatically  mapped,  as  a  function  of  rigidity. 

This  cutoff  rigidity  mapping  technique  utilizes  the  trajectory 
parameterization  approach  to  locate  regimes  of  rigidity  within  which 
trajectory  approach  is  forbidden  due  to  sane  particular  planet  (or 
atmosphere)  intersection.  By  examining  the  details  of  the  trajectory  at 
the  point  where  it  makes  the  intersection  it  is  possible  to  efficiently 
establish  the  rigidity  width  of  the  resulting  penunbral  band,  whether  it  be 
fine  or  wide.  By  automating  strategies  by  which  this  examination  is  made, 
the  program  very  rapidly  determines  the  position  and  width  of  the  penunbral 
bands,  obviating  the  process  of  tracing  trajectories  at  very  fine  rigidity 
intervals  that  previously  had  been  necessary  to  map  the  penumbra  in  any 
direction.  In  the  new  technique,  when  a  wide  penunbral  band  is  responsible 
for  the  failure  of  trajectories,  it  is  necessary  only  to  calculate  a  few 
trajectories,  widely  spaced  in  rigidity,  in  order  to  determine  the  band 
position  and  width.  On  the  other  hand,  narrev  bands  are  reliably  detected 
and  their  details  established  to  an  equivalent  precision. 


7.2  POSSIBLE  USE  OF  MAPPED  FTNUMBRAL  DETAIL  TO  ESTABLISH  WHETHER  OOEMIC 
PAY  OBSERVATIONS  COULD  BE  USED  TO  CHECK  AND  3MFR0VE  THE  CURRENT 
MATHEMATICAL  MODELS  OF  THE  GEOMAGNETIC  FIELD 
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Seme  of  the  penumbral  features  explored  exhibited  a  very  sensitive 
dependence  on  the  geomagnetic  field  model.  For  example,  the  penumbral 
"islands"  discussed  in  Chapter  4  shew  such  a  sensitive  dependence.  The 
possibility  of  using  these  structures  to  experimentally  determine  the 
appropriateness  of  any  given  mathematical  model  of  the  Earth's  field,  for 
example,  exists.  In  principal  if  a  detector  could  be  vised  to  observe 
primary  cosmic  rays  arriving  at  a  site  in  the  direction  of  the  axis  of 
symmetry  about  which  the  islands  are  disposed,  the  spectral  characteristics 
of  the  observed  flux  ought  to  reflect  the  characteristics  of  the  islands  - 
the  rigidity  at  which,  for  example,  the  penumbral  island  finally  disappears 
with  reducing  rigidity  ought  to  be  visible  as  a  discontinuity  in  the 
spectrum  of  the  cosmic  ray  flux.  Because  this  rigidity  value  is  very 
directly  dependent  on  the  field  modeling,  a  comparison  oculd  be  made 
between  the  cosmic  ray  flux  in  the  actual  situation  and  that  predicted. 

Figure  7.1  shows  how  the  "islands"  and  the  surrounding  major  penumbral 
structure  change  form  with  rigidity.  The  penumbral  transmission  in  the 
axial  direction  of  the  islands  would  have  the  form,  as  a  function  of 
rigidity,  shown  in  figure  7.2.  Such  a  spectral  signature  should  be 
detectable.  Another  phenomenon  which  oould  be  observable,  too,  in 
principle,  is  the  gradient  associated  with  the  structure  in  the  north-west, 
which  is  a  constant  position  over  a  range  of  rigidity. 

A  major  difficulty  arises,  however,  in  that  it  would  be  necessary  to 
locate  the  detector  in  fixed  locations  at  relatively  lew  satellite 
altitudes  in  order  to  accumulate  the  required  net  detector  observations. 

The  alternative  to  this  impossible  requirement  would  be  to  carry  cut  a 
very  long  term  experiment,  and  to  correlate  the  events  observed  at 
different  locations  with  the  penumbral  structure  at  the  site  where  the 
event  was  registered,  and  by  shifting  the  relative  rigidity  scales  and 
effectively  superimposing  events,  accumulate  a  sufficiently  great 
normalized  event  total  to  check  the  form  of  the  calculated  spectra,  and 
thus  the  field  model. 

Ground  level  observations  are  unsuitable  for  two  reasons.  First,  the 
flux  within  the  atmosphere  would  consist  of  secondary  particles,  and  so  the 
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FIGURE  7.1 


Variation  of  the  size  and  shape  of  the  penuntoral  islands  and 
surrounding  major  penumbral  structure  with  changing  rigidity.  The  rigidity 
values  range  from  7.1  to  7.9  Gv  in  o.l  GV  steps.  The  parameters  relating 
to  these  plots  are  the  same  as  for  figure  5.5  (10°  N  latitude,  270°  E 
longitude,  400  km  altitude;  IGRF80  field) .  These  are  zenith-azimuth  plots 
of  the  same  form  as  figure  5.5,  and  may  usefully  be  ocxipared  with  that 
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rigidity  spectrum  of  the  primaries  would  not  be  directly  observable.  In 
any  case,  the  penumbral  islands  lie  at  high  zenith  angles,  and  would  be 
below  the  horizon  at  ground  level. 


Penumbral  transmission  in  the  axial  direction  (28°  zenith,  235° 
azimuth)  of  the  penumbral  islands  in  figure  7.1,  as  a  function  of  rigidity 
(hatching  represents  forbidden  trajectories) . 


An  alternative  is  to  use  an  experiment  on  board  a  geosynchronous 
satellite  to  correlate  the  existing  cosmic  ray  fluxes  with  those  predicted 
using  particular  field  models.  In  order  to  pursue  this  possibility 
further,  the  parameterization  technique  was  used  to  map  the  penumbral 
structure  as  it  would  exist  at  the  orbit  of  a  stationary  satellite.  Figure 
7.3  presents  the  penumbral  maps.  These,  calculated  for  a  range  of  rigidity 
values,  shew  essentially  the  shadow  of  the  Earth  at  the  given  rigidities. 

A  great  degree  of  sensitivity  does  exist  in  seme  of  these  structures.  In 
particular,  at  the  rigidity  where  the  high  rigidity,  essentially  line-of- 
sight  pattern  bifurcates  into  two  forbidden  regions  (corresponding  to  the 
Earth's  shadow  as  seen  by  particles  approaching  from  both  directions  along 
the  local  field  line  bundle)  the  pattern  displays  great  sensitivity  to 
change  in  the  calculational  parameters.  An  experiment  on  board  a  satellite 
in  a  stationary  orbit  could  be  used  to  monitor  the  cosmic  ray  fluxes  in 
this  direction,  and  to  correlate  the  directional  intensities  with  field 
model  and  time  (a  time  dependence  due  to  the  magnetospheric  compression 
would  be  expected) . 
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FIGURE  7.3 


structures  for  a  range  of  rigidity  tmluas,  •Mff'cptate 
s  geosyndrerou.  satsllit.  lying  .t  Utitui.  0°,  lcrgitud.  0°. 
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